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Abstract 

A review of the author's results is given. Inversion formulas and stability re- 
sults for the solutions to 3D inverse scattering problems with fixed energy data 
are obtained. Inversion of exact and noisy data is considered. The inverse po- 
tential scattering problem with fixed-energy scattering data is discussed in detail, 
inversion formulas for the exact and for noisy data are derived, error estimates 
for the inversion formulas are obtained. The inverse obstacle scattering problem 
is considered for non-smooth obstacles. Stability estimates are derived for inverse 
obstacle scattering problem in the class of smooth obstacles. Global estimates for 
the scttering amplitude are given when the potential grows to infinity in a bounded 
domain. Inverse geophysical scattering problem is discussed briefly. An algorithm 
for constructing the Dirichlet-to-Neumann map from the scattering amplitude and 
vice versa is obtained. An analytical example of non-uniqueness of the solution to a 
3D inverse problem of geophysics and a uniqueness theorem for an inverse problem 
for parabolic equations are given. 
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1 Introduction 



In this paper 3D inversion scattering problems with fixed-energy data are discussed. 
These problems include inverse problems of potential, obstacle, and geophysical scattering 
(IPS, IOS, IGS). 

Inverse potential scattering problem is discussed in detail: uniqueness of its solution, 
reconstruction formulas for inversion of the exact data and for inversion of noisy data 
are given and error estimates for these formulas are obtained. These estimates yield the 
stability estimates for the solution of the inverse scattering problem. 

For the inverse obstacle scattering the uniqueness theorem is proved for rough do- 
mains, stability estimates are obtained for C 2,x domains, < A < 1, that is, for domains 
whose boundary in local coordinates is a graph of C 2,x function. Reconstruction formulas 
are discussed. 

For inverse geophysical scattering the inverse scattering problem is reduced to inverse 
scattering problem for a potential. 

Construction of the Dirichlet-to -Neumann map from the scattering data and vice 
versa is given. Analytical example of nonuniqueness of the solution of an inverse 3D 
problem of geophysics is given. 



The results discussed in this paper were obtained mostly by the author, see [l|], [10]- 



however the presentation and some of the estimates are improved in this paper. 
Only some selected results from the cited papers are included in this review. 

1.1 The direct potential scattering problem. 

We want to study the inverse potential scattering problem of finding q(x) given some 
scattering data. 

Consider the direct scattering problem first and let us formulate some basic results 
which we need. 
Let 

[V 2 + k 2 - q{x)]u{x, a, k) = in M 3 , x G M 3 , (1.1) 



u = e lkax + A(a',a,k)— + o(-\ r := \x\ -> oo, a' :=- (1.2) 

r \r J r 

Here u(x, k) is the scattering solution, k = const > is fixed. Without loss of 
generality we take k — 1 in what follows unless other choice is suggested explicitly. A 
unit vector a G S 2 is given, where S 2 is the unit sphere in M. 3 . Vector a has a physical 
meaning of the direction of the incident plane wave, while a' G S 2 is the direction of the 
scattered wave, k 2 is the fixed energy. The function A(a', a, k) is called the scattering 
amplitude. It describes the first term of the asymptotics of the scattered field as r — > oo 
along the direction a' — -. 
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The function q(x) is called the potential. We assume that 
qeQ :=Q a nL°°(M 3 ), 



Qa ■= {q ■ q{x) = q(x), q(x) G L (B a ), q(x) = if \x\ > a}, 



;i.3) 



where a > is an arbitrary large fixed number which we call the range of q(x), and the 
overbar stands for complex conjugate. 

In many results q G Q a is sufficient, but q G Q is used in the proof of a crucical 
estimate (2.17) below. 

1.2 Review of the known results. 

Let us formulate some of the known results about the solution to problem (1.1)-(1.2), 
the scattering solution. These results can be found in many books, for example, in the 



appendix to JT0|, where a brief but self-contained presentation of the scattering theory is 
given. 

1.2.1 The scattering problem has a unique solution if q G Q a . 

In fact, the above result is proved for much larger class of q ([[J, |J), but for inverse 
scattering problem with noisy data it is necessay to assume q(x) compactly supported 
IPH - Indeed, represent the potential q(x) as q = q\ + q<i, where q\ — for \x\ > a and 



qi = q for \x\ < a. Call q2 the tail of the potential q. If one assumes a priori that 
q = 0(\x\~ b ), where b > 3, then the contribution of the tail of the potential to the 
scattering amplitude is of order 0(|a| 3_fc ) and tends to as a — > oo. At some value of a, 
say at a = a , this contribution becomes of the order of the noise in the scattering data. 
One cannot, in principle, discriminate between the noise and the contribution of the tail 
of the potential for a > clq. Therefore the tail of q for a > ao cannot be determined from 
noisy data. 
One has 

sup \u(x,k)\ < c, k = const > 0. (1.4) 

aiGK 3 

By c > we denote various constants. If q G Q a then u(x, k) extends as a meromorphic 
function to the whole complex fc-plane. Let G(x,y, k) denote the resolvent kernel of the 
self-adjoint Schrodinger operator Lu = —V 2 + q(x) in L 2 (R 3 ): 

(L — k 2 )G(x, y, k) = 5(x — y) in M 3 , (1.5) 



lim 



o\x\ 



ds = 0, y is fixed , k > 0. (1.6) 
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The function u(x, k) can be defined by the formula: 

e ik\y\ / i \ y 

G(x,y,k) = —r-u(x,a,k) + o — , — = -a, (1.7) 

4?r M \\y\J \y\ 

where o = O {^2 j is uniform with respect to x varying in compact sets and formula 



'1.7) can be differentiated with respect to x [1J, [16 



The function G(x, y, k) is a meromorphic function of k on the whole complex fc-plane. 
It has at most finitely many simple poles ikj, kj > 0, 1 < j < J in C + := {k : Irak > 0} 
and if q(x) ^ 0, q G Q a , infinitely many poles, possibly not simple, in C_ = C \ C+. 
There are no poles on the real line except, possibly at k — 0. 

The functions Uj(x) G L 2 (M 3 ), solving (1.1) with k = ikj, are called eigenfunctions 
of the discrete spectrum of L, —kj are the negative eigenvalues of L. There are at most 
finitely many of these if q G Q a . 

The eigenfunction expansion formulas are known: 

/(,■) V. /•»,(./•) • /' f(Z)u(x,Z)d£, \£\ = k, £ = ka, 

where 

1 



fj : = (/,«i)L 2 (K 3 ), f(0 = TTTTs / f{x)u(x,a,k)dx, 
(see e.g. |TJ). 

If is the resolution of the identity of the selfadjoint operator L, and E\(x, y) is its 
kernel, then 

— A [ - ^ = —ImG(x, y, y/X) = — — g / u(x, a, VX)u(y, a, VX)da, A > 0. 
dX 7r \Qtx a J S 2 

1.2.2 Properties of the scattering amplitude 



The scattering amplitude has the following well-known properties (see e.g. [10 
A(a', a, k) = A(— a, —a' , k) (reciprocity), 



A(a', a, k) — A(a', a, —k), k > (reality), 
k 



ImA(a',a,k) = — / A(a' , (3, k)A(a, (3, k)d/3, k>0 (unitarity). 
4tt J S 2 

In particular, 

k f 

ImA(a,a,k) = — / \A(a, (3, k)\ 2 df3 (optical theorem). 
4tt J S 2 
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If q G Q a , and k — 1, then the scattering amplitude is an analytic function of a' and a 
on the algebraic variety 

3 

M := {0 : <E C 3 , 0-9 = 1}, -w:=J2 e j w r ( L8 ) 

This variety is non-compact, intersects M 3 over S* 2 , and, given any ( G R 3 , there exist 
(many) 0,0' E M such that 

0'-0 = f, |0|->oo, 0,0' EM. (1.9) 

In particular, if one chooses the coordinate system in which £ = te 3 , £ > 0, e3 is the 
unit vector along the a^-axis, then the vectors 

0' = ^e 3 + C 2 e 2 + Ciei, = ~e 3 + C 2 e 2 + Ciei, Ci + C 2 2 = 1 - j, (1-10) 

satisfy (1.9) for any complex numbers Ci and £ 2 satisfying the last equation in (1.10) and 
such that |Ci| 2 + IC2I 2 — ► 00. There are infinitely many such (1X2 E C If q G Q a than 
the function A(a', a, k) is a meromorphic function of G C which has poles at the same 
points as G(x, y, k). 
One has 

-AnA(a', a, k) = / e~ iha '' x q(x)u(x, a, k)dx. (1.11) 
The ^-matrix is defined by the formula 

S = I + —A, S*S = SS* = I, (1.12) 

27T 

and is a unitary operator in L 2 (S 2 ). 
If Vg G Q a then 

0(x, a, fc) := e~ lka ' x u(x, a, k) = 1 + — — / 9(3; - ra)dr + o f \ J , fc — > 00. (1.13) 

22K Jo V^/ 



Therefore 



q(x) =a-V x lim {2ifc[0(x, a, k) - 1]} , (1.14) 



and 

4 frv. rv. fc") = — 

4:71 



A(a,a,k) = —— I q(x)dx + o(l), fc — > +00. (1-15) 
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1.2.3 The fundamental equation. 

Denote u + := u(x, a, k),u~ := u(x, —a, —k), k > 0. Then u + = Su~ , that is 



ik f 

u + = u~ H / A(a',a,k)u~(x,a',k)da'. (1-16) 

2tt J S 2 

1.2.4 Completeness properties of the scattering solutions. 

a) Ifh(a) G L 2 {S 2 ) and 

/ h(a)u(x, a, k)da = Vi G B' R := {x : \x\ > R}, k > is fixed (1-17) 
t/ien /i(a) = 0. 

Let N D (L) = {w : Lw = in D, w G iJ 2 (D)}, where D C f 3 is a bounded domain, 
H 2 (D) is the Sobolev space. 

b) The set {u(x, a, k)}y a( - s2 is total in N D (L — k 2 ), that is, for any e > 0, however 
small, and any fixed w G N D (L — k 2 ), there exists v e (ot) G L 2 (S 2 ) such that 

— / u(x, a, k)v £ {a)da^H 2 (D) < £■ (1-18) 
Js 2 ' 

The z/ e (o;) depends on iy(x). 

1.2.5 Special solutions 

There exists ^(x, 9,k), <E N D (L — k 2 ), such that 

[V 2 + A; 2 - g(x)]^ = in R 3 , ijj = e ik6 ' x [l + R(x, 6, k)], 6 G M, (1.19) 

fin 101)3 

\\R\\l°°(d) < c K 1 ' , |0| -> oo, G M, (1.20) 
1 6> 1 2 

II^IU^d) < , |0|-oo, 9 <E M, (1.21) 
where Z) C M 3 is an arbitrary bounded domain. 
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1.2.6 Property C for the pair {L\ — k 2 ,L 2 — k 2 } 

Let Ljw := Yl\nt\=o a jm(x)d^w(x), x E W 1 , n > 2, j = 1,2, be linear formal partial 
differential operators, that is, formal differential expressions. 

Let Nj = Nju(Lj) := {w : LjW = in D}, where D C R n is an arbitrary fixed 
bounded domain and the equation is understood in the sense of the distribution theory. 
Consider the subsets of Nj, j = 1, 2, which form an algebra in the sense that the products 
W\W 2 E LP (D), where Wj E Nj, p' = and 1 < p < 00. If p = 1 define p' = 00, and if 
p = 00 define p' = 1. We write Via,- E Nj meaning that Wj run through the above subsets 
of A^, 

Definition 1.1. We say that the pair of linear partial differential operators {Li,L 2 } has 
property C p if and only if the set {wiu^} is total in L P (D), that is, if f(x) E L P (D) and 



/ f(x)w 1 (x)w 2 (x)dx = Vwj E Nj,j = 1, 2, (1.22) 

then 

f(x)=0. (1.23) 
// the above holds for any p > 1, we say that property C holds for the pair {L\, L 2 }. 

Theorem 1.1. 1.1 / fli|/ . Let Lj = —V 2 + Qj{v), qj{x) E Q a , k = const > is arbitrary 
fixed. Then the pair {Li — k 2 , L 2 — k 2 } has property C. 

Proof. Note that ipj E Nj, j = 1, 2, where ipj are defined in section 1.2.5 above. Without 
loss of generality take k = 1, let ip(x, 9, 1) := ip(x, 9). One has 

Vi(x, 9')^ 2 {x, -9) = e ^'-^(l + R X ){1 + R 2 ). 

Choose 9' ,9 E M such that (1.9) holds with an arbitrary fixed ^6l 3 . Then 

V>i(z, 9')ip 2 (x, -9) = e^ x {\ + o(l)) as \9\ 00. (1.24) 

Since the set {e^' x }v£gR3 is total in L P (D), p > 1, D C E n is a bounded domain, the 
conclusion of Theorem 1.1 follows. □ 



Remark 1.1. One cannot take unbounded domain D in the above argument because o(l) 
in (1.24) holds for bounded domains. 

One can take the space of f(x) larger than L l (D), for example, the space of distribu- 
tion of finite order of singulartiy if q(x) is sufficiently smooth [IT] . 

Theorem 1.2. The set {ui(x, a, k)u 2 (x, (3, fc)}v«,/?e5 2 ,fc>o complete in L P (D), where 
D C IR 3 is an arbitrary fixed bounded domain, and p > 1 is fixed. 

Proof. The conclusion of Theorem 1.2 follows from Theorem 1.1 and (1.18). □ 
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1.2.7 Properties of the Fourier coefficients of A(a',a). 
We denote A(a! , a, k)\k=i '■= A(a',a), and write 

oo „ 

A(a',a) = ^2A e (a)Yt(a'), A e (a) :=/ A(a',a)Yi(a')da', (1.25) 

where Yt(a') = Y^ m (a'),—£ < m < £, summation over m, is understood in (1.25) and 
similar formulas below, e.g. (1.31), (1.37), etc, 



Yt_ m (a) 



(2£ + l)(£-m)\ 
(£ + m)\ 



e^P, m (costf), 



Y e , m {a) = (-l) l+m Y^ m ( a ), Y ltm {-a) = (-l)% m (a). 

Here P^^cosi?) = (sin ^) m d ^ e Q ^^ , < m < £, Pe(x) is the Legendre polynomial, (i?, <p) 

are the angles corresponding to the point a G S 2 , P^_ m (cos $) = (— l) m ^~^| P^ m (cost?), 
< /// < /. 

Consider a subset M' C M consisting of the vectors = (sin$cos</9, sin d sin ip, cost?) 
where d and ip run through the whole complex plane. Clearly G M, but M' is a proper 
subset of M. Indeed, any 9 G M with 63 7^ ±1 is an element of M' . If 03 = ±1, then 
cos i? = ±1, so sin-0 = and one gets 9 = (0,0, ±1) G M' . However, there are vectors 
9 = (61,62,1) G M which do not belong to M' . Such vectors one obtains choosing 
0i,02 G C such that 6\ + 6\ = 0. There are infinitely many such vectors. The same is 
true for vectors (6±, 62, —1). Note that in (1.9) one can replace M by M' for any £ G M 3 , 
^2e 3 . 

Let us state two estimates (|fL5ll): 

™v\M<*)\<c(^Qf\ (1.27) 



and 



1 „r\Im9\ 

\Y e (6)\< Vr>0, 6eM', (1.28) 

V47T |M r )l 



where 

Mr) ■= (^f J ^ {r) = 271 1 (I)' 11 + aS ^ °°' (L29) 

and J^(r) is the Bessel function regular at r = 0. Note that Yg(a'), defined by (1.26), 
admits a natural analytic continuation from S 2 to M by taking d and p in (1.26) to be 
arbitrary complex numbers. The resulting 6' G M' C M. 
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1.2.8 A global perturbation formula. 

Let Aj(a',a) be the scattering amplitude corresponding to qj G Q a , j = 1,2. Define 
A := A x - A 2 , p := gi(ar) - q 2 {x). Then Q 

— 4tcA(cx', a) = / p(x)iti(x, a)u2(x, —a')dx. (1.30) 

1.2.9 Formula for the scattering solution outside the support of the poten- 
tial. 

Let supp (q) C B a . The fixed-energy scattering data A(a', a)Va', a, or, equivalently, the 
data {A^(a)}^ = o,i,2,...Va G S" 2 , allow one to write an analytic formula for the scattering 
solution u(x, a) in the region B' a := M. 3 \B a : 



e i a -x + J2A i (a)Y i (a')h i (r), r := \x\ > a, a' := -, (1.31) 



where Ae(a) are defined in (1.25), Yf(a') are defined in (1.26), 

Mr) :=e'H'«)^/|;<,M, 

H^\r) is the Hankel function, and the normalizing factor is chosen so that 

h e (r) = — [l + o(l)l asr^oo. (1.32) 
r 

Note that [2, formula (7.1463]: 

\Hf\r)\ 2 = — / K (2rsht)(e m + e~ m )dt, (1.33) 

JO 



where sht := e *~ 6 . This formula implies that |^(r)| is a monotonically increasing 
function of t. 

It is known [Q], formula 8.478] that r|/i£(r)| 2 is a monotonically decreasing function 
of r if £ > 0, and 

i e \/2 f2(V 

h e {r) = — - [l + o(l)], i^+oo, r>0. (1.34) 

r \er J 

The following known estimate can be useful: 

\Y t , m {a)\ <c£^-\aeS 2 , (1.35) 
where Y^ m (a) are the normalized in L 2 (S 2 ) spherical harmonics (1.26). 
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Let us give a formula for the Green function G(x,y, k) (see (1.5), (1.6)) in the region 
\x\ > a, \y\ > a, where suppg(a;) C B a . Let g(x,y,k) := f~nrz^i an d denote by Aye the 
Fourier coefficients of the scattering amplitude: 



A(a',a) = J2M«)YtW) = J2 Ae e Y e (a)Y e (a'). 



;i.36) 



1=0 



e',e=o 



Then 





G(x,y,k) = g(x,y,k) + — ^ Ae> e Y e/ (-a)Y e (a')he(k\x\)ht>>(k\y\), \x\ > a, \y\ > a, 



e\e=o 



(1.37) 



where a' := A, a := A. 

fI' Ivl 

Indeed, clearly the function (1.37) solves (1.5) in the region \x\ > a, \y\ > a, where 
q(x) =0, it satisfies (1.6), and 



G(x,y,k) 



,ik\y\ 



Air \y | 



e i ka - x + k 2 A e , e Y e ,(a)Y e (a')h 



tWx\ 



£',e=o 



+ o 



1.38) 



as \y\ — > cxd, — = —a. 



By (1.36), (1.31), (1.25), and (1.7), it follows that the function (1.37) has the same 
main term of asymptotics (1.38) as the Green function of the Schrodinger operator. 
Therefore the function (1.37) is identical to the Green function (1.5)-(1.6) in the region 
\x\ > a, \y\ > a. 

2 Inverse potential scattering problem with fixed- 
energy data 

The IPS problem can now be formulated: given A(a', a) Vet', a G S 2 , find q(x) G Q a . 
Throughout this section k — 1. 



2.1 Uniqueness theorem. 

The first result is the uniqueness theorem of Ramm [|13| , |14 



Theorem 2.1. Ifqi,qz G Q a and Ai(a', a) = A 2 (a',a) Ma' G Sf, Wa G S%, where S 
j = 1,2, are arbitrary small open subsets of S 2 , then q±(x) = (j^O 3 ?)- 
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Proof. The function A(a',a) is analytic with respect to a' and a on the variety (1.8). 
Therefore its values on Sf x S 2 extend uniquely by analyticity to M x M. In particular 
A(a',a) is uniquely determined in S 2 x S 2 . By (1.30) one gets: 

/ p(x)ui(x, a)U2(x, — a')dx — Va, at' G S 2 . (2.1) 

JB a 

By property C (section 1.2.6), formulas (1.22) -(1.23)) and by (1.18), the orthogonality 
relation (2.1) implies p(x) = 0. □ 

2.2 Reconstruction formula for exact data. 

Fix an arbitrary (el 3 and choose arbitrary 9', 9 satisfying (1.9). 
Denote 

q(0 ■= [ e-* x q(x)dx. (2.2) 

JB a 

Multiply (1.11) by u(a, 9) G L 2 (S 2 ), where is(a, 9) will be fixed later, and integrate 
over S 2 with respect to a: 

-47r / A(a',a)v(a,9)da = / e " 4a '' x / u(x,a)u(a,9)daq(x)dx. (2.3) 
Js 2 JB a Js 2 

If q G Q a , then estimates (1.27) and (1.28) imply that the series (1.25) converges, 
when a' is replaced by 9' G M, uniformly and absolutely on S 2 x M c . where M c C M 
is an arbitrary compact subset of M. Formula (1.11) implies that a' can be replaced by 
9' G M, since B a is a compact set in M 3 . 

Define 

p(x) := p(x; v) := e~ t8 " x / m(x, a)i/(a, - 1, (2.4) 

Js 2 



and rewrite (2.3), with a' = 9', as 



-Air A(9',a)v(a,6)da= t~ lB - x+l °- x [p(x) + l]q(x)dx 

JS 2 JB a 

— q{Q + / e~ % ^" x p{x)q(x)dx = q + e, 

JB a 



(2.5) 



where 

|e| < ||g|| ||p||o, h\\a ■= \\q\\v*(B a )- (2-6) 

The following estimate (see [|T5|], estimate (2.17) and its proof in section 6 below) 
holds for a suitable choice of u(a, 9): 

a<c\9\~\ |0|->oo, 9eM. (2.7) 
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From (2.5) and (2.7) one gets the reconstruction formula for inversion of exact fixed- 
energy 3D scattering data: 



lim 1-471-/ A(9',a)v(a,9)da\ (2. 



0,0'eM. 

and the error estimate: 

-4tt J A(9\a)v(a,9)da = q(£) + o(J-^ , \9\ oo, 9 E M, (2.9) 

where (1.9) is always assumed. 

Let us give an algorithm for computing the function u(a,9) for which (2.7), and 
therefore (2.9), hold, given the scattering data A(a',a) W, a G S* 2 . 

Fix arbitrarily two numbers and 6 such that 

a < a x < b, (2.10) 
and define the L 2 -norm in the annulus: 

llpf := / \p\ 2 dx. (2.11) 

J ai<\x\<b 

Consider the minimization problem 

M:=d(9), (2.12) 



where the infimum is taken over all v G L 2 (S 2 ). 

It is proved in JOd] (see also section (6.3) below) that 

d{9) < c|0| _1 if 9 G M, \9\ > 1. (2.13) 

The symbol \9\ ^> 1 means that \9\ is sufficiently large. The constant c > in 
(2.13) depends on the norm ||g|| a but not on the potential q(x) itself. An algorithm for 
computing a function u(a,9), which can be used for inversion of the fixed-energy 3D 
scattering data by formula (2.9), is as follows: 

a) Find any approximate solution to (2.12) in the sense 

\\p(x,u)\\ <2d(9), (2.14) 

where in place of 2 in (2.14) one could put any fixed constant greater than 1. 

b) Any such z/(a, 9) generates an estimate of g(£) with the error o (j§\^)i 1^1 > °°- 
This estimate is calculated by the formula 

g:=_47r/ A{9',a)u(a,9)da, (2.15) 
Js 2 

where u(a,9) G L 2 (S 2 ) is any function satisfying (2.14). 
We have obtained the following result: 
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Theorem 2.2. One has 

\n - nl.f )\ < ^ 



q-m\<^, |0|-oo, 9 EM, (2.16) 



provided that (2.14) an d (1-9) hold. 

Proof. The proof is the same as the proof of (2.5) - (2.7) and is based on the following 
estimate 0, 111 : 



U < c + I6")- 1 ) , |0|>1, 6 E M. (2.17) 
The proof of (2.17) is not simple ||15|| . It is given in section 6. □ 



2.3 Stability estimate for inversion of the exact data. 

Let the potentials q E Q, j — 1, 2, generate the scattering amplitudes Aj(a / , a). 
Let us assume that 

sup |Ai(a',a) - A 2 (a',a)\ < 5. (2.18) 

a',a€S 2 

We want to estimate p(x) := qi(x) — q2{x). 
The main tool is formula (1.30). 
The result is ([0, 



Theorem 2.3. If qj E Q and (2.18) holds th 



en 



sup |&(0 - g 2 (£)| < c^fl, 6 -> 0, (2.19) 

where £o > is an arbitrary large fixed number and the constant c > does not depend 
on 5 > 0, 5 — > 0. 

Proo/. Multiply both sides of (1.30) by ^(a, 0)i/ 2 (-a', 2 ), where 0j e M, j = 1,2, 
0i + 02 = £, |0i| — ► oo, and integrate with respect to a and a' over S" 2 , to get: 

— An / v4(a', a)^i(a, 0i)^ 2 (— a', 2 )dada' = 

Js2 Js2 (2.20) 

dxp(x) / «i(a;, a)^i(a, 0i)c?a / u 2 (x, /3)z^ 2 (/3, &2)dP, (3 = —a 1 . 
B a Js 2 Js 2 



Choose v\ and z/ 2 such that 

111 < 



WpM)\\ < TTtV ^-eM, (2.21) 



j 1 
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where Pj{vj) '■= pj(x,Uj) := e ldjX J g2 Uj(x, a)uj(a, 9j)da — 1, and note that i|4 — > 1 as 
|#i| — > oo, ^i,^2 G 6>i + 6> 2 = £, |£| < £o, £o is an arbitrary large but fixed number. 
From (2.18), (2.20) and (2.21) one gets 

< c (l^r 1 + c^HU^IN^)) , (2.22) 



where 9 := 9\. One can choose v\ and v% such that ([0, see also section 6.7 below) 

||^(a,^-)|U W ^ce^lM*!, = 1; |^| _ oo, 9 2 = £ - 9, (2.23) 

where c > stands for various different constants. 
Thus (2.22) yields: 

sup |p(£)| < cmin [s" 1 + Cl 5e C2Sln|s| l , < 5 « 1, (2.24) 

Kl<& s>>1 

where c, ci and C2 are some positive constants, s := \9\, 8 <C 1 means that 5 > is small 
and s 3> 1 means that s > is large. However, our argument is valid for s > 1 and 
0<5<±. 



One gets 



min Is- 1 + Cl 5e C2Sln|s| l := V {5) < c 3 ^|^, 5-0, (2.25) 
s>o L J mo 



and the minimizer is 



S = S (5) = c 2 - 1 i i^L[l + o(l)], 5-0. (2.26) 
From (2.24)-(2.26) one gets (2.19). □ 



Remark 2.1. In the above proof the difficult part is the proof of (2.23). Estimate (2.23) 
can be derived for v(a, 9) such that 

r e~~ bK 

\\i[)(x,9)- u(x,a)u(a,9)da\\ L 2 (Bb) < c , « := \Im9\, 9 £ M, 9 > 1, (2.27) 

Js 2 K 

and \\v(a,9)\\i,2(s 2 ) — inf ; where the infimum is taken over all v £ L 2 (S 2 ). 

In section 6.7 we consider the problem of finding v £ L 2 (S 2 ) with minimal norm 
||z/|| L 2( 5 2) := a(u) among all u(a,9) which satisfy the inequality: 

\\ip(x) — / u{x, a)v(ot, 9)doc^L 2 (B b ) < £■ (2.28) 
Js 2 
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The necessity to consider the v with the minimal norm Hz/H^^) comes from the 
simple observation: there exists a sequence of v n e L 2 (S 2 ), ||i/ n || l 2 (s 2 ) — 1; such that 

u(x,a)u n (a)da\\u2/ Bj) \ — > as n — > oo. (2.29) 

s 2 

To prove (2.29) note that 

a) = - JL——^ q (y)u(y, a)dy := e ia ' x - Tu, (2.30) 

where the operator (J + T) -1 := I + T% is a continuous bijection of C(B b ) onto itself, and 
C(Bb) is the usual space of continuous in i??,, b > a, functions equipped with the sup- 
norm JT5|. Since T is compact in C(Bb), the above statement follows from the injectivity 
of I + T, which we now prove: 

If / + Tf = 0, then / is extended to C(IR 3 ) by the formula / = — Tf, and satisfies 
the following equation (V 2 + l — q(x))f = in M 3 and the radiation condition of the type 
(1.6) with k = 1. Therefore f(x) = and the injectivity of / + T is proved. 

Thus I + T and J + T\ are continuous bijections of C[B^) into itself for any b > a. 

Writing u(x, a) = (I + Ti)e ia ' x , one concludes that (2.29) is equivalent to 

|| / e m ' :r z/ n (a)(ia|| L 2( jBf) ) — > as n — > oo. (2.31) 
is 2 

Existence of a normalized sequence z/ n (a) satifying (2.31) follows from the compactness 
of the operator 

Q : L 2 (S 2 ) -> L 2 (5 6 ), := / e 4a '^(a)da. 

is 2 

Of course, the same argument is applicable to the operator Qiv := f s2 u(x,a)b , (a)da, 
but the bijectivity of I + T in C(Bb), b > a, is of independent interest. 

It follows from (2.29) that, for a given e > 0, one can find z/ in (2.29) with an arbitrary 
large norm ||z/|| i 2( > 52- ) . By this reason we are interested in v with minimal norm. Estimate 
(2.23) gives a bound on the growth of the minimal value of the norm ||v|| £,2(152), where 
v = i/(a,d) satisfies (2.28) with e = « := \Im9\, \6\ —>■ 00, 6 E M. 



2.4 Stability estimate for inversion of noisy data. 

Assume now that the scattering data are given with some error: a function A$(a',a) is 
given such that 

sup \A{a',a) - A s {a',a)\ < 5. (2.32) 

a' ,aGS 2 

We emphasize that Ag(a',a) is not necessarily a scattering amplitude corresponding 
to some potential, it is an arbitrary function in L°°(S 2 x S 2 ) satisfying (2.32). It is 
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assumed that the unknown function A(a', a) is the scattering amplitude corresponding 
to a q E Q. 

The problem is: Give an algorithm for calculating q such that 



sup \q~q(0\ <V(&), ?7(<5)^0 as 5^0, 
Kl<& 



(2.33) 



where £ > is an arbitrary large fixed number, and estimate the rate at which r)(5) tends 
to zero. 

An algorithm for inversion of noisy data will now be described. 
Let 



N(5) :-- 

where [x] is the integer nearest to x > 0, 

N(S) 



\\n5\ 
ln|ln5| 



(2.34) 



A s (0',a):=f]A se (a)Y e (d'), A se (a) := I A s (a',a)Y e (a')da', (2.35) 

7^o Js2 



N(S) 



u s (x, a) := e ia - x + ^ A 5e (a)Y e (a')h e (r), 



[2.36) 



ps(x; v) := e l9 ' x / u s (x, a)v(a)da - 1, 9 e M, 
Js 2 



(2.37) 



^(5) : = e-T^W, 7 = const > 0, 

a(u) := || z/|| £2(s2), k := |/m#|. 
Consider the variational problem with constraints: 

|0| = sup := &(8), 

\9\ [\\ps(v)\\+a(is)e Kb fi(5)] < c, 9 e M, 
where the norm is defined in (2.11), and it is assumed that 

0'-e = £, 9,9' G M, 



^2.38) 
^2.39) 

^2.40) 
(2.41) 

^2.42) 
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where £ G M 3 is an arbitrary fixed vector, c > is a sufficiently large constant, and the 
supremum is taken over 9 G M and v G L 2 (S 2 ) under the constraints (2.41). 

Given (gK 3 one can always find 9 and 9' such that (2.42) holds. 

We prove that ~&(S) — > oo, in fact 



&(6) > c- 



hi6\ 



Let 9(5), h , s(a) be any approximate solution to problem (2.40)-(2.41) in the sense that 



|0OS)| > 



d(5) 



Calculate 



Qs 



—Att / As(9', a)vs(oi) da. 

I S 2 



(2.44) 



(2.45) 



Theorem 2.4. One has 



sup \q s -q(0\ < c 
l€l<& 



(ln|ln<5|) ; 
|ln5| 



as <5 — > 0. 



Proof. 



& - q(0 = - ^ A s (9', a)u s (a)da - q(£) 
Js 2 



-An f A(9', 
Js 2 



a)vs(oi)da — q(£) 



+ 4tt / [A(9',a) - A 5 (9',a)]u 5 (a)da. 
Js 2 



The rest of the proof consists of the following steps: 
Step 1. We prove that 



\\p(x,is 5 )\\ < c[\\ Ps (x,u s )\\ +a(v 5 )e Kb fi(5)] < 



where the norm is defined in (2.11). 

This estimate and (2.43) imply (see the proof of (2.19) and (2.26)) that 



-471-/ A(9', a)v s (a)doL - q(£) 
Js 2 

Step 2. We prove that 



< c 



(ln|ln5|) : 
\\n5\ 



l S 2 



A(9' , a) - A s (9' , a u s (a)da 



, , h , , c (In I ln5|) ; 
< ca(v 5 )e Kb n(5) < < c- V 



ln5| 



(2.46) 



(2.47) 



(2.48) 



(2.49) 



(2.50) 
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where 9' = 6'(8) = £ + 9(5), and the pair {9(5), vs(a)} solves (2.40)-(2.41) approximately 
in the sense specified above. (See formula (2.44)). 

This estimate follows from (2.41) and from the inequality 



\\A(9',a) - A 5 (9',a)\\ LHs2) < ce Kb fi(5). 
Let us prove (2.51). One has 

N(5) 

\\As(9'(5),a)-A(9'(5),a)\\ LHs2) < || £ [A sl (a) - A(«)] Y e (9')\\ LHs2) 



(2.51) 



£=0 



(2.52) 



+ || A e (a)Y e (9')\\ LHs2) :=I 1 + I 2 , 

£=N(S)+l 

where N(5) is given in (2.34) and & := 9'(5). Using (1.28), (1.29) one gets 



h < c5N 2 (5) 



e Kb (2N) N ^ +1 
(ehjW) ' 



(2.53) 



Here we have used the estimate 



sup ^2 \ Ae(oc) - Ase(a) 



< 4tt5 2 , 



which follows from (2.32) and the Parseval equality and implies 



sup 

aeS 2 



A s (a) - A 5l (a) 



(2.54) 



We also took into account that there are (N + l) 2 spherical harmonics Y e = Y^ m with 

< t < N, because YL=-t 1 = 2^ + 1, and Eio( 2 ^ + l ) = ( N + l f- For lar S e N one 
has (N + l) 2 = iV 2 [l + o(l)], N -> oo, so we write (N + l) 2 < cN 2 , c> 1. 
To estimate I 2 , use (1.27) -(1.29) and get: 

JV+1 



V27/ 



< ce 



Kai 



Minimizing with respect to N > 1 the function 



: (2N) 



N+l 



5N 2 K ^\{„ +s N+ \ 0<s :=-<!. 



(eft) 



A? 



a 
ai 



one gets 



mm 

N>1 



5N 



2 (2N) N+1 



(eb) 



N 



+ S 



N+l 



< ;;e ~7iV(<5) = c/i (5) ; 7 = i n n > 0, 



ai 



(2.55) 



(2.56) 



(2.57) 



where iV(<5) is given in (2.34) and fi(5) is defined by (2.38). Thus, from (2.52) -(2.56) one 
gets (2.51). Theorem 2.4 is proved. □ 
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2.5 Stability estimate for the scattering solutions. 

Let us assume (2.18) and derive the following estimate: 
Theorem 2.5. // q\,q 2 G Q a and (2.18) holds then 

sup \u\{x, a) — u 2 (x, a) \ < c/i(<5), \x\ > a, (2.58) 

x£B' a 

aes 2 

where fi(5) is defined by (2.38) and (2.34), an d c > is a constant. 
Proof. Using (1.31), one gets: 

oo 

u x - u 2 = ^[A a {a) ~ A e2 (a)]Y e (a')h e (r), r > a. (2.59) 

£=0 

As stated below formula (1.33), one has 

\he(r)\ < \h e+1 (r)\ , r > 0. (2.60) 

From (2.58) one gets: 

oo „ b 

\\ui-u 2 \\ 2 L 2 {Bb \ Bai) = ^2\A ei (a)-A i2 (a)\ 2 r 2 \h e (r)\dr, a<a x <b. (2.61) 

It follows from (2.60) and (1.34) that 



(2.62) 



rb rb 

sup / r 2 \he(r)\ 2 dr = / r 2 \h N (r)\ 2 dr < 
0<£<N J ai J a x 

f2N\ 2N f b dr (2N\ 2N _ 2N 

From (2.60)-(2.61) one gets 

K " Mh(B b \B a ) < ^ (^f) + c E (l^il 2 + M , (2.63) 

where we have used the monotone decrease of r\h^{r)\ 2 as a function of r. Using estimate 
(1.27)in order to estimate |A^(o;)|, j = 1,2, one gets: 

/2N\ 2N ( a \ 2N 
hi - u 2 \\ 2 L2{Bb ^ Bai) < c5 2 1^— j +c[^-J , «<«i- (2-64) 
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Minimization of the right-hand side of (2.63) with respect to iV > 1 yields, as in (2.56), 
the estimate similar to (2.56): 

hi-^2\\l HBb \ Bai) <ce-^ 5 \ 7 = ln^>0. (2.65) 

Since u\ — u<i := w solves the equation 

(V 2 + l)w = in B' a := R 3 \B a , (2.66) 

one can use the known elliptic estimate: 

IMIff2(Di) < c (||(V 2 + l)w\\ L 2 (D2) + \\w\\ L 2 {D2) ) , Di C D 2 , (2.67) 

where D\ is a strictly inner subset of D2 and c = c(Di, D2), and get: 

Mh^d,) < ce~^ 5 \ (2.68) 

where Di is any annulus a\ < a 2 < \x\ < a 3 < b. By the embedding theorem, (2.68) 
implies (2.58) in R 3 . □ 

2.6 Spherically symmetric potentials. 

If q(x) = q(r), r := \x\, then 

A(a',a) = A{a' ■ a), A e {a) = A e Y e {a). (2.69) 



In [^3| the converse is proved: if q E Q a and (2.69) holds then q(x) = q(r). 

The scattering data A(a', a) in the case of the spherically summetric potential is 
equivalent to the set of the phase shifts 5g. The phase shifts are defined as follows: 

1 + —A e = e 2i5e , A t = 4vre^ sin^, k = 1. (2.70) 

From Theorem 2.1 it follows that is q = q(r) E Q a then the set {^}^=o,i,2,... determines 
uniquley q{r). A much stronger result is proved by the author in ||24|| . To formulate this 
result, denote by £ any subset of nonnegative integers such that 



E 7 = o °- ( 2 - 71 ) 



Theorem 2.6. (WM) Ifl{. x ) — ?( r ) ^ Qa then the data {^}v£eL determine q(r) uniquely. 

In |] an example is given of two quite different potentials qi(r) and q2(r), piecewise- 
constant, qj(r) = for r > 5, for which 

5^ - 5 p) 



sup 

£>0 



< 10 and sup \qi(r) — q2(r)\ > 1, where q\ and g 2 

r>0 



are of order of magnitude of 1. 

This result shows that the stability estimate (2.19) is accurate. 
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3 The direct obstacle scattering problem with fixed- 
frequency data 

Let D C I 3 be a bounded domain (an obstacle) with boundary S. We want to study 
the obstacle scattering problem under minimal smoothness assumptions on S. 

Recall that if S is C 1,A , < A < 1, that is, in local coordinates S is a graph of a 
C l,x - function x% = g(x'), x' := (xx,X2), g G C l,x , then the obstacle scattering problem 
consists of finding the scattering solution u(x, a), which satisfies the equations: 

(V 2 + l)u = in D' := R 3 \D, (3.1) 
Tu = 0, (3.2) 

u = e ia ' x + A(a'a) — + o ( - ) , r := \x\ -»• oo, - = a', (3.3) 
r \r J r 

where a G S 2 is given and (3.2) is the Dirichlet condition if Tu = u, the Neumann 
condition if Tu = un, or the Robin condition if Tu := un + cr(s)u, a(s) G L°°(S), 
Ima(s) = 0, iV is the exterior normal to S. We took the wavenumber k = 1 without loss 
of generality. 

If S is very rough (non-smooth), N may be not defined on S. The minimal assump- 
tions on the smoothness of S, under which the existence and uniqueness of the solution 
to the direct scattering problem are established, were introduced in jl7[] and [TI|]. These 
assumptions are: 

Ax) If Tu = u then D is an arbitrary bounded domain, that is a bounded open set. 
A 2 ) If Tu = un then the assumption on S is: 

i : H l (D' a ) — > L 2 (D' a ) is compact , D' a := D' D B a . (3.4) 

Here i is the embedding operator, i/ 1 is the Sobolev space, a > is such a number that 
the ball B a := {x : \x\ < a} contains D, and D' := M. 3 \D. 
A 3 ) If Tu = un + cr(s)w, then the assumption on S is: 

% : H\D a ,) -> L 2 (D a i) and : ^^D^) -> L 2 (S) are compact. (3.5) 

Here the integration on in the definition of L 2 (S) is understood with respect to the 
two-dimensional Hausdorff measure on S. 

The usual classes of domains in the theory of Sobolev spaces are: 

1) domains satisfying the cone condition, 

2) Lipschitz domains, 
and 

3) extension domains, (see [0, |8|]). 
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They all are such that the above assumptions A 2 ) and A3) hold. Let us recall the 
definitions of these domains: D satisfies the cone condition if each point of D is the 
vertex of a cone contained in D along with its closure, the cone in the local coordinates 
is the region x' 2 < ax 2 , < X3 < b, and a, b > are fixed positive constants. A domain 
D is Lipschitz if each point of 5* := dD has a neighborhood [/cl n , such that U H D 
can be mapped onto a cube by a quasi-isometric map. 

A homeomorphic / : D x —>■ D 2 is called quasi-isometric if 

r i/(s)-/(s )i LT^W^MI <- ^ 

hm sup : : < M, hm sup ; ; < M 

x^x \x — yo\ y^yo \y — yo\ 

for any x G D\ and any y G D 2 and the Jacobian det f'(x) preserves its sign in D\. 

A domain D is an extension domain in H 1 if there exists a linear continuous operator 
E : H\D) -> ^(M 3 ), Eu = u on D for all u G H l (D). 

An extension domain may fail to satisfy the cone condition, but a bounded domain 
satisfying the cone condition is the union of a finite number of extension domains. 

We will use also the domains with finite perimeter. A domain D has finite perimeter 
if and only if || Vxd(s ; )||_bv(ir 3 ) < °°j where Xd(x) is the characteristic function of D and 
BV is the space of functions of bounded variation . 

Let us state the result from [|18|1 . 

Theorem 3.1. Problem (3.1) - (3.3) has a unique weak solution if: 

a) D is an arbitrary bounded domain (open set) in R 3 and Tu = u, 

b) D is a bounded domain, condition A 2 ) holds and Tu = un 

c) D is a bounded domain, condition A3) holds and Tu = Un + o~(s)u, a(s) G L°°(S), 
Ima(s) = 0. 

The solution to (3.1) - (3.3) for rough boundaries is understood in the weak sense and 
has the following properties: 

if Tu = u then u EH 1 (D' R ) n C% C {D'), u G L 2 (d 1 , j^s), a > 1, and (3.3) holds 



a 

with 



A(a',a) = -— I e- ia '- s u N (s,a)ds, (3.6) 

47T 



s 



where H 1 (D' R ) is the Sobolev space of functions which vanish on S, 

b) if Tu = u N , then u G H l (D' R ) n C^ C (D'), u G L 2 (d', y^), a > 1, and (3.3) 
holds with 

1 f de' ia ' s 

A{a,a) = — I u(s,a)ds, (3.7) 

J s 
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c) if Tu = u N + a(s)u, then u G H\D' R ) n C^ C (D'), u e L 2 (p' , i+j^pr), a > 1, and 
(3.3) holds with 



y4(c/, a) = — 



47T 



<9e 



9iV, 



+ cr(s)e" 



u(s, a)ds. 



(3.8) 



In JTjJ and JTB| Theorem 3.1 was proved with the operator V 2 replaced by a general 
second-order selfadjoint elliptic operator. The weak solution is defined in the case Tu = 
Un + cr(s)u, as a function u G H 1 ^^) fl C^ C (D') for any R > Rq, B Ra D D, which 
satisfies (3.3) and satifies the integral relation: 



(VuV0 - u(f))dx - / a(s)ucj)ds = V0 G Hq(D' 

D> Js 



(3.9) 



Here Hq(D') is the set of if 1 (D') functions vanishing near infinity and cis is the two- 
dimensional Hausdorff measure on S. Formula (3.9) makes sense for domains D with 
finite perimeter. 



3.1 Uniqueness theorem for inverse obstacle scattering. 

The inverse obstacle scattering problem consists of finding S and the boundary condition 
on S given A(af, a) W G S 2 . The scattering amplitude A(a', a) in the obstacle scattering 
problem satisfies conditions listed in section 1.2.2. In particular, it admits unique analytic 
continuation from S* 2 x S* 2 onto M x M. 

Let us outline the proof of the uniqueness theorem for inverse obstacle scattering 
problem with fixed-frequency data. This theorem belongs to the author [ 16| , but we give 
a new proof |2tJ , see also Jl9| . 



Theorem 3.2. If A(a', a) = ^2(0/, a) Va ; , a, running through arbitrary small open sub- 
sets of S 2 , then D\ = D2 '■— D, and the boundary condition on S := dD are uniquely 
determined. 

Proof. As in the proof of Theorem 2.1, the data determine uniquely Aj(a', a) on S 2 x S 2 
so that Ai(a', a) = A 2 (a', a) Va', a G S 2 . If one has already proved that D\ = D 2) that 
is, 5*i = 5*2 := S, then the boundary condition on S is uniquely determined because the 
scattering solution u(x, a) is uniquely determined by the scattering amplitude in D' (and 
is analytically determined by formula (1.31) in B a /). Thus, the limiting values of — on S 
are uniquely determined. If the limit is zero (almost everywhere on S) then Tu = un, if 
the limit is infinity (almost everywhere on S) then Tu = u, and if the limit is a function 
—a(s), then Tu = + a(s)u. Therefore the main point is to prove that S is uniquely 
determined by A (a', a) Va', a & S 2 . Let us prove this. 

Assume the contrary: St ^ S 2 . Let D 12 := D 1 U D 2 , S 12 = dD 12 , D 12 := D x n D 2 . 
Denote by D\ a connected component of D\ 2 \D 2 . We want to show that D\ 2 \D 2 is an 
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empty set. An important tool is the formula [T2I| similar to (1.30): 
- Air \Ai(a\ a) - A 2 (a', a)] = 

[miat(s, — a')u 2 (s } a) — wi(s, — cx')u 2 n(s } «)] ds. 



(3.10) 



Si 2 



This formula holds for domains with finite perimeter. 
If A x = A 2 W, a G S 2 , then (3.10) yields: 



/ [um(s, — a')u 2 (s, a) — U\(s, — o!)v^n(s, &)] ds = Vo/, a E S 2 . (3-H) 
From (3.11) and (1.7) one derives: 

/ [G 1N {x, S )G 2 (s,y)-G 1 {x,s)G 2N ( S ,y)}ds = Vx,y e D' 12 . (3.12) 

JSl2 

From (3.12) and Green's formula one gets: 

G l (x,y) = G 2 (x,y) Vx, y G D' 12 . (3.13) 
This leads to a contradictionjrnless D\ = D 2 . 

Indeed, if T)\ ^ D 2 , i.e., if D\ is not empty, take a point s on the boundary of T)\ 
which belongs to S\. This point is an interior point for D' 2 . Thus 

G 2 (s, y) — >• oo as y — > s. (3-14) 

On the other hand, since s G Si one has TGi(s,y) = 0, i.e., Gi(s,y) = if Tu = u, 
GiN(s,y) = if Tu = un, Giw(s,y) + a(s)Gi(s,y) = if Tu = un + cr(s)u. In all the 
three cases 

TG 2 (s, y) — >• oo, as y — > s, (3.15) 

while 

r(?i(s,y) =0. (3.16) 
From (3.13), (3.15) and (3.16) one gets a contradiction. Theorem 3.1 is proved. □ 



Remark 3.1. In fftdy a uniqueness theorem is proved for inverse obstacle problem with 
the transmission boundary condition: the boundary condition and the boundary and the 
wavenumber in the interior of the obstacle are uniquely determined by the fixed-frequency 
scattering data A(a',a) Vcn', a G S 2 . 

Remark 3.2. The basic reason to consider the domains with finite perimeter is the va- 
lidity of the Green formula in such domains j^j. 

In it is proved that is Di and D 2 are domains with finite perimeter then Di is 
also such a domain. 
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3.2 Stability estimate for inverse obstacle scattering. 

Consider the scattering problem (3.1) - (3.3) with Tu = u, for example. Let D\ and 
D 2 be two arbitrary obstacles in the class O x which consists of the bounded domains 
whose boundaries can be covered by finitely many balls Bj, on the patch Sj = SOBj the 
boundary is described in the local coordinates by the equation x 3 = gj(x'), x' = (xi, x 2 ), 
where gj(x') e C 2 > x , < A < 1, and 

sup ll^llco^fs,) < c , (3.17) 
3 

where cq does not depend on the choice of D e O x . 

Let Aj(a', a) be the scattering amplitude corresponding to the obstacle Dj, j = 1, 2. 

Assume that (2.18) holds. Define the symmetric Hausdorff distance between Si = dDi 
and 5*2 = dD 2 : 

p := p{5) := max{sup inf \x — y\, sup inf \x — y\}. (3.18) 

xeS 2 ^ eS i xeSi v^ s 2 



The basic stability estimate [2^| can now be formulated. Let A m (a', a) be the acattering 
amplitude corresponding to D m and Tu = u on S m := dD m , m = 1,2. 

Theorem 3.3. // (2.18) and (3.17) hold then 

P(5)<c,(^Y, 5^0, (3.19) 



ln<5| 

where ci,c 2 are positive constants independent of5>0. 

Proof. Let us sketch the steps of the proof. 
Step 1. 

p(5) -»• as 5 -> 0. (3.20) 

This follows from the uniquness theorem 3.2 and from the compactness of the set O x in 
the space C 2,fJl , p < A. For simplicity of the presentation we assume that j = 1, that is, 
there is just one patch in the covering of S, for example, S is star-shaped. 
Step 2. There exists an integer m such that 

cd m (x) < \w(x)\ < ce cdCix \ w(x) := u x {x, a) - u 2 (x, a), (3.21) 

where d(x) is the distance from a point x G D' l2 to Si 2 , we assume that d(x) ~ p, 
e = ce~ lN ( s \ 7 = const > 0, N{5) = lr [ , and c > here and below stand for various 
constants indepent of 6 and x, and we assume that d(x) ~ p. Symbol d(x) ~ p means 
that c±d(x) < p < c 2 d(x) with some constants ci,c 2 > 0. 
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Let us show that (3.21) implies (3.19). From (3.21) one gets, replacing d by p, 
dropping w and taking log: 



hip < c + cp c In e. 



(3.22) 



Recall that c stands for different constants. 
From (3.22) one gets 



ln± 

p 



< 



ln± 

£ 



(3.23) 



Since p' 3 < as p 



and (5 > 0, estimate (3.23) implies 

'ln|hi5| XC2 



p < c 



ln± 



Cl 



In 5 1 



where we have used the definition of e, namely e = ce ~< N ( S ) ^ which implies -^r ~ \\^s\ • 
Let us give the details of the proof. 

Step 1. Assume that (3.20) is false. Then p{5) > c > as 5 — ► 0. Since O x is a 
compact set in C 2,/i , p < A, one can select sequences Si™^) and S , 2n (5) which converge in 
C 2 ' M to Si and £2 correspondingly, as 5 — > 0. Since A(a', a) depends continuously on S 
(see 



27 , ilbl) in the sense 



lim sup \A n (a', a) — A(a', a)\ = 0, 



(3.24) 



where the limit is taken in the process S n — > 5 in C 2,M , < p < A, one concludes that 
Ai(a', a) = A 2 (a', a) for the limiting surfaces Si and S2. By the uniqueness theorem 3.1 
it follows that Si = 5*2. Therefore p := p(Si, S2) = 0. However, p = lim^^o p(8) > c > 0. 
This is a contradiction which proves (3.20). 

Step 2. The function w(x) := Ui(x, a) — m 2 (x, a), where Uj(x, a) is the scattering 
solution corresponding to the obstacle Dj, j = 1, 2, solves the equation 



(V 2 + l)w = in D[ 2 , Di 2 := D x U D 2 , 



satifies the radiation condition 



lim 

r— >oo 



dw 



d\x\ 



iw 



ds = 0, 



(3.25) 



(3.26) 



and 



\w\\ C 2,x {D > l2 ) < C. 



(3.27) 
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It is proved in 0, vol.3, p. 14, that solutions to elliptic second order equations with smooth 
coefficients cannot have zeros of infinite order up to the boundary without vanishing 
identically. This implies existence of an integer m > for which the left inequality (3.21) 
holds. 

The proof of the right inequality requires some preparations. 
Let us sketch the steps of this proof. 

Step 3. By the result of Theorem 2.5 one gets estimate (2.58) in the region \x\ > a, 

Although in Theorem 2.5 the functions u\ and u 2 were the solutions to the Schrodinger 
equations with potentials vanishing in B' a , the estimate (2.58) is proved for any solutions 
to equation (3.25) whose difference satisfy (3.26) and (2.18). In particular, estimate 
(2.58) holds for our w. Let us define e = cfi(6) = ce^ N{ - 5 \ N{8) = Jj^L. 

Step 4. Let us prove the right inequality (3.21). Extend w from D' l2 into D 12 so that 
the estimate similar to (3.27) holds: 

IMIc 2 .*(R3) < ci. (3.28) 

This is possible: using, for example, the known Stein's theorem one can extend Uj into 
Dj since Sj is C 2,x - smooth, j = 1,2, and then w := U\ — u 2 will be the C 2,x - smooth 
extension of w from D' 12 into Di 2 . Define 

f( x ) ■= -(V 2 + l)w in R 3 , (3.29) 

where we denoted by w the extended to R 3 function w(x). Then f(x) = in D' 12 , 
f(x) E C A (R 3 ), and 

w(x) = [ - J (y)dy, x e R 3 . (3.30) 

Jd 12 4:Tr\x-y\ 

Denote \x\ = r, \y\ = p. Set z := re %Lp . Choose any point x G D[ 2 , in a neighborhood 
of Di, a connected component of D 12 \D 2 , such that d(x) := dist(x, Si) ~ p := p(S). By 
step 1, p(5) — > as 5 — > 0, so that p(8) is small for small 5. 

Consider analytic continuation ofw(x), defined by (3.30), on the complex z-plane (as 
in [11]). Let u be the angle between vectors x and y, x G D' 12 . Since 

\x — y\ = y/(r — pe lu) ){r — pe~ iU} ), 

this analytic continuation, that is, replacement r — > z, in the expression \x — y\, is possible 
if re l{p — pe luJ ^ and re lLp — pe~ luJ ^ 0. Choose a point O on Si closest to x and the 
coordinate system in which the origin is at O, and the xix 2 — plane is tangent to Si 
at the point O. Since Si is sufficiently smooth, there exists a cone K with an opening 
6*o > and vertex at O which belongs to D' l2 and its axis passes through the point x. 
The function \x — y\ admits analytic continuation in the z-plane from the ray r > to 
the sector \cp\ < 9q. Since there are no points of Di 2 inside the cone K, the expression 

( ■= [(re ilp - pe iu )(re iip - pe~ iu )]^ , Im( > 0, r > 0, | arg^| = |(^| < 6> , 
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does not vanish in the region r > 0, | axgz\ := \(p\ < 6 . Therefore the function (3.30), 
which we denote W(z), considered as a function of z, admits analytic continuation in the 
sector | argyj| < 9q, r > 0, and satisfies the following inequalities there: 

\W{z)\ < c, | argz| < O , r > 0, (3.31) 
|W(z)|<e, arg^ = 0, r > a, (3.32) 

where 



One can map the sector | argz| < 9 , conformally onto the half-plane, z — > t = z™o , 
| argt| < |. Then W(z) := v(t), where v(t) is analytic in the half-plane | argt| < |, and 
satisfies there the inequalities 

\v(t)\<c, |argt|<~, |v(t)|<e, t > a x . (3.34) 

The known two-constants- theorem (see 0) and (3.34) imply: 

\v(t)\< Cl e h ®, (3.35) 

where h(t) is the harmonic measure corresponding to the domain R on the complex 
plane t with the boundary consisting of the lines {Ret = 0, — oo < Imt < oo} and 
{Imt = 0, Ret > d\ > 0}. Recall that h = h{t\, t 2 ), the harmonic measure, is a harmonic 
function which solves the problem: 

Ah := h tltl + h t2t2 = in R, (3.36) 

h = at t± = 0, — oo < £2 < 00, h = 1 at t 2 = 0, t\ > a, 

h(ti,t 2 ) is bounded at infinity, t — t\ + %t 2 . 

By the maximum principle, 1 > h(ti,t 2 ) > in R, and, by the Hopf lemma, (see || 
p.34), 

dh 

— > at tt = 0. (3.37) 

ut\ 

Thus 

h(t u 0) > cti, < ti < t , c = const > 0, (3.38) 

where to > is a sufficiently small number. 

From (3.35) and (3.38) it follows that in a sufficiently small neighborhood of the origin 
one has |t>(ti,0)| < ce ctl . Returning to the ^-variable one gets 

\W(z)\ < ce c|2|c , z = r=\x\>0, (3.39) 

where c > stands for various constants. 

Since W(r) = w(x), inequality (3.39) is identical to the right inequaltiy (3.21). 
Theorem 3.2 is proved. □ 
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Remark 3.3. An interesting open problem is to construct S in the inverse obstacle scat- 
tering problem analytically from noisy data Ag(a',a) in the way it is done in section 
(2.4), formula (2.45), for the reconstruction of the potential in the inverse potential scat- 
tering problem, or even from exact data in the way it was done in section 2.2, formula 
(2.8). In the next section we prove that such a reconstruction formula does exist for exact 
data for the inverse obstacle scattering problem. 



3.3 Existence of a reconstruction formula. 

Let A(a', a) be the scattering amplitude corresponding to an obstacle D, and assume, 
for example, that Tu = u, so that the Dirichlet condition holds on S. Take a' = 9' G M 
in (3.6), multiply (3.6) by a v(a,9) and integrate with respect to a over S 2 to get: 

-4tt / A(9',a)u(a,9)da = e~ ie ' s u N (s,a)u(a,9)dads. (3.40) 
Js 2 Js Js 2 

Here u(a, 9) G L 2 (S 2 ) is some function which is chosen so that (3.41) holds. We prove 
below the following lemma: 

Lemma 3.1. The set {un(s, a)}\/ a€ s 2 is total in L 2 (S). 

This implies that, given an arbitrary small number rj > and 9 G M, there exists a 
i/„(a,0) G L 2 (S 2 ) such that 

f de w ' s 
|| J ^u N (s,a)u v (a,9)da- -^-\\ L 2 (S) <r]. (3.41) 

From (3.40), with v = v v , and (3.41) one gets: 

-47rlim / A(9',a)v r) (a,0)da= / e^'' 5 — — ds (3.42) 
v^oJ S 2 ' Js oN 

Let us assume 9,9' G M, 9 — 9' = £, £ G M 3 is an arbitrary fixed vector. 
By Green's formula one gets: 

■ \ L - -f L - -f-«>. <-) 

where Xd(x) — 1 in D, Xd(x) = in D'. 

From (3.42) and (3.43) one gets an inversion formula for finding D from the data 
A(a', a): 

Xd(0 = 87r|£r 2 lim / A(9', a)v n {a, 9)da. (3.44) 
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Before proving the basic Lemma 3.1, a remark is in order: in contrast to our theory for 
inverse potential scattering problem (see the inversion formula (2.8) and the algorithm 
for calculating the function v(a, 9) in (2.15)) we do not give an algorithm for calculating 
the function u v (a,9) in (3.44). 

Finding such an algorithm is an open problem. 

We now prove Lemma 3.1: assume the contrary. Then there exists a function / G 
L 2 (S) such that 

f{s)u N {s,a)ds = Va 6 S 2 . (3.45) 

Js 

Define 

v(y):= [ f(s)G N (s,y)ds, (3.46) 

where G(x, y) is the Green function: 

(V 2 + l)G(x, y) = —5(x - y) in D', (3.47) 



G{x,y) = 0, xeS, (3.48) 



G satisfies the radiation condition (1.6). 
We claim that (3.45) implies 



v{y) = in D', (3.49) 



where v(y) is defined in (3.46). 
Indeed, 



(V 2 + l)v = in D', (3.50) 

and 

v(y) = o ( — J as \y\ — > oo. (3.51) 

\\y\J 

The relation (3.51) follows from (3.45) and (1.7). 

It is well known pf, p. 25, that (3.50) and (3.51) imply (3.49). 

If (3.49) holds, then, taking y — > a G S (along the normal to S at the point a) in 
(3.46), one gets 

f(s) = 0. (3.52) 

Indeed, if o G S then 

G N (s, y) -> 5 s (s - a) as y -> a, (3.53) 
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where 5s(s — a) is the delta-function on the surface S (see formula (3.58) below). 
From (3.53), and (3.49) and (3.46) formula (3.52) follows. 
Let us give a proof of (3.53). 
Consider the problem: 

(V 2 + l)v = in D', (3.54) 
v = f(s) on S, (3.55) 

v satisfies the radiation condition (3.26). (3.56) 
This problem has a unique solution representable by the Green formula: 

v(x) = / f(s)Gjsr(s,x)ds. (3.57) 
Js 

Since (3.55) holds for this unique solution (3.57) taking x — > a 6 S one gets (3.53), 
where 5s{s — a) is the distribution which acts by the formula 

/ f(s)5 s (s-a)ds = f(a). (3.58) 

This argument requires f(s) to be continuous if one understands the delta-function 
in the usual sense. However, if Ss(s — a) is understood as the kernel (in the distributional 
sense) of the identity operator in some space of functions f(s), for which problem (3.54)- 
(3.56) has a unique solution, then (3.58) makes sense in the space for example, in L 2 (S), 
and (3.58) is understood in this case as equality of the elements of this space. In the 
case of L 2 (S) this means that the equality holds almost everywhere on S with respect to 
two-dimensional Hausdorff measure on S. Lemma 3.1 is proved. 

Thus, formula (3.44) is proved. □ 

4 Limiting procedure and stability estimates 

Let D C 1R 3 be a bounded domain, x{ x ) be the characteristic function of D, t > be 
a parameter, q(x) = t\{x) be the potential. Consider the potential scattering problem 
(1.1)— (1.2). We prove that the scattering solution u(x, a;t) converges, as t — > +oo, to 
u(x, a), the scattering solution corresponding to the obstacle D, and u(x, a) satisfies the 



Dirichlet boundary condition on S := dD. This result is old ||29|| . We also prove the 



following, more recent estimates [21 



c c 
\u(x } a;t)\\ L 2 {D) < —, \\u(x } a;t) -u(x,a)\\ L2(3l) < (4.1) 
£2 £2 
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||V«(a;,a;f)|| £2(Bil) < c(R), B R D D, (4.2) 

Q 

\\u(x,a;t)\\ L 2 {s) < — , (4.3) 

where D' is a strictly inner compact subdomain of D' := ~§L 3 \D. 

Assume there are two obstacles Di and D 2 , Dj C 2 ' x , < A < 1. Let 

A(a', a; ti, t 2 ) :— A 1 (a', a; ti) - A 2 (a', a; t 2 ). 

Then we prove the following stability estimate: 

sup \A(a',a;t u t 2 )\ < c I + p J , (4.4) 

where t± 2 = min(ti,t 2 ) and p is the symmetric Hausdorff distance between Si and 5*2, 
defined by (3.18). 

Let us prove the above estimates. 

One has 

[V 2 + 1 - tx(x)] u(x, a; t) = in M 3 . (4.5) 

Define 

ll9 f \u(x)\ 2 dx _ . * 

u : = / 7TT .2^ ' (7>1 - ( 46 ) 
Jrs (1 + |a;| 2 ) 2 

We drop the a-dependence in a; t). 
Assume that 

\\u(x;t)\\ <c, (4.7) 

where c = const > does not depend on t. 

As t — > +oo, one can select, using (4.7), a weakly convergent in the norm (4.6) 
sequence, denoted again u. Thus 

M(x,t)^M. (4.8) 

By elliptic estimates (see formula (6.1) below), formulas (4.8) and (4.5) imply 

\\u(x;t)\\ H2 ^ < c, \\u(x,t) - u(x)\\ H 20 l) -> as t -> oo. (4.9) 
Multiply (4.5) by u(x,t) and integrate over a ball -Br D £>, to get 

/ [|Vn(x; t)| 2 + t)| 2 ] drr = f \u(x; t)\ 2 dx + / J^y«ds (4.10) 
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where the bar stands for complex conjugate. 
From (4.10), (4.7) and (4.9) one gets 



[ \Vu(x;t)\ 2 dx <c(R), [ \u(x;t)\ 2 dt<^-. (4.11) 
Jb r Jb r t 

This yields (4.2) and the first inequality in (4.1). 

Let us prove (4.3). The embedding theorem yields (see, e.g., [[7], p. 66): 



\u 



(x;t)\\ L 2 (s) <c(e\\Vu{x;t)\\ L 2 {D) +e- 1 \\u{x;t)\\ L 2 {D) ) , < e < 1. (4.12) 



Take e = t 5 and use (4.11). Then (4.12) yields (4.3). Let us prove the second inequality 
(4.1). 

Denote v := u(x; t) — u(x). Then 

V 2 v + v = in D'; v satisfies (3.26), (4.13) 

v = u(x;t) on S. (4.14) 



By Green's formula, one gets 



v(x)= [ u(s;t)^^-ds, (4.15) 
Js 9N S 

where G(x,y) is the Green function for the Dirichlet operator V 2 + 1 in D' . 

\v{x)\ < \\u(s;t)\\ L 2 (s) \\ dG ^ X) \\ L 2 is) . (4.16) 

From (4.16) and (4.3) the second estimate (4.1) follows. 

It is clear that the function u = u(x) in (4.8) is the scattering solution u(x, a) corre- 
sponding to the obstacle scattering problem with the Dirichlet condition on S. Indeed, 
u = on S (see(4.3)), u solves equation (4.13) in D', u — e ia ' x satisfies the radiation 
condition. These three conditions determine u uniquely, and since u(x, a) satisfies these 
conditions, it follows that u(x,a). 

Let us prove (4.4). 

If D\ = D2 = D, then formula (1.30) yields 

\A(a',a;t 1 ,t 2 )\ < \t 2 - h\ Wu^x, a; ti)\\ L 2( D) \\u 2 (x, -a'; t 2 )\\^(D) < c-^ — (4.17) 

(tih) 2 

If t\ — t 2 — +00, then formula (3.4) yields: 

\A{a',a)\ ^cpfaSt). (4.18) 
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In the general case a combination of the above estimates yields (4.4) (see PT| ) an d also 
the estimate 



\A(a', a; £ x , t 2 )\ < c [\h - £ 2 | + p(S u S 2 )} 



(4.19) 



useful when \t\ — t%\ is small and £i,£2 £ [l,£o], where £0 > 1 is fixed. 
If Dt = D 2 , then (4.19) yields 



\A(a', a;£i,£ 2 )| < c\t\ - t 2 



(4.20) 



5 Inverse geophysical scattering with fixed-frequency 
data 

Consider the problem 



where v(x) G L 2 (JR 3 ) is a compactly supported real- valued function with support in R_, 
the lower half-space. In acoustics u has the physical meaning of the pressure, v(x) is the 
inhomogeneity in the velocity profile. We took the fixed wavenumber k — 1 without loss 
of generality. The source y is on the plane P := {x : X3 = 0}, i.e., on the surface of the 
Earth, the receiver 16P. 

The data are the values {w(x, y)}\/ x ,yeP- 

The inverse geophysical scattering problem is: given the above data, find v(x). 

The uniqueness theorem for the solution to this problem is obtained in [pOj] , |i~l| . 

Problem (5.1) -(5.2) differs from the inverse potential scattering by the source: it is 
a point source in (5.1) and a plane wave in (1.2). Let us show how to reduce the inverse 
geophysical scattering problem to inverse potential scattering problem using the "lifting" 
|nj , p5| . Suppose the data w(x,y), x G P, y G P, are given. Fix y and solve the 
problem: 




(5.1) 




(5.2) 



(V 2 + l)w = in R 3 



{x : x 3 > 0} 



(5.3) 



w = w(x,y), 



x G P, 



(5.4) 



w satifies (5.2). 



(5.5) 
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This problem has a unique solution and there is a Poisson-type analytical formula for 
the solution to (5.3)-(5.5), since the Green function of the Dirichlet operator V 2 + 1 in 
the half-space is known explicitly, analytically: 

e i\x—y\ „ik\x—y\ 

Gi(x, y) = — 1 - — — , y := (yi,y 2 , -ys)- (5.6) 

Att\x — y\ Atx\x — y\ 

Therefore the data w(x,y)\/x G P determine uniquely and explicitly (analytically) 
the data w(x,y)Wx G Ri, y G P. We have lifted the data from P to as far as 
^-dependence is concerend and get w(x,y)Wx,y G given w(x,y)\/x,y G P. 

If w(x,y) is known for all x, y G M+, then one uses formula (1.7) and calculates 
analytically the scattering solution u(x, a) corresponding to the potential q(x) := —v(x) 
and k — 1, where a G S* 2 := {a : a G S 2 ,a^ < 0}. Given u(x, a) for all x G 
and a G S* 2 , one can calculate the scattering amplitude A(a',a) W G S* 2 := {a : a G 
S 2 ,a 3 >0}. 

If the scattering amplitude A(a', a), corresponding to the compactly supported q(x) = 
—v(x) G L 2 (M 3 ) is known W G 5^., Vet G S* 2 , then the uniqueness of the solution to 
inverse geophysical problem follows from theorem 2.1. 

Stability estimates obtained for the solution to inverse potential scattering problem 
with fixed-energy data remain valid for the inverse geophysical problem: via the lifting 
process one gets the scattering amplitude A(a', a) corresponding to the potential q(x) = 
—v(x), and the stability estimates for q(£), obtained in sections 2.3 and 2.4, yield stability 
estimates for v(£) = —q(£). 

Practically, however, there are two points to have in mind. The first point is: if the 
noisy data u$(x,y) are given, where sup x yeP \us(x , y) — u(x,y)\ < 5, then one has to 
overcome the following difficuly in the lifting process: data ip(x,y),x,y G P, such that 
\ip(x,y)\ < 5, may not decay as \x\ — ■> oo, \y\ — > oo on P, and this brings the main 
difficulty 

The second point is: if one uses the inversion algorithms presented in sections 2.2-2.4, 
then one uses the data A(a', a)\fa! ', a G S* 2 . Of course, the exact data A(a', a)Wa' G 
S* 2 , a G S 2 , determine uniquely the data A(a', a) Wa', a G S* 2 , but practically finding the 
full data from the partial data is an ill-posed problem. 

6 Proofs of some estimates. 

Here we prove some technical results used above: estimates (1.18), (1.19), (1.20), (1.30), 
(2.13), and (2.17). 

6.1 Proof of (1.18). 

It is sufficient to prove (1.18) with L 2 (D) in place of H 2 (D): since w(x) and 

u(x, a)v e {a)da 



J 

Js 2 
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solve equation (1.1), the elliptic estimate (see [Q) 

IMItf 2 (.Di) < c [||L^|| L 2 (I52) + ||<^|| Z 2 (D2) ] , Di C 5 2 , (6.1) 

where L = V 2 + 1 — q{x), Di is strictly inner subdomain of D 2 and c = c(D 1 ,D 2 ) = 
const > 0, implies that IMIi^rDi) < c IMU 2 0D 2 ) if Lip = 0. If (1.18), with k = 1 and 
L 2 (D) in place of H 2 (D), is false then 



w(x) u(x : a)u(a)da = \/u(a) G L (S ). (6.2) 



Therefore 



w(x)u(x, a)dx = Va G S 12 . (6.3) 



This implies 



/ w(x)G(x, y)dx = \/y e D', 
Jd 



(6.4) 



where G(x, y) is the Green function of the operator L. 

Indeed, denote the integral on the left-hand side of (6.4) by (p(y). Then 

Lip = in D'; ip = o ( — ] as \y\ — > oo. (6.5) 

The second relation (6.5) follows from (6.3) and (1.7). From (6.5) one gets (6.4) by 
lemma 1 on p. 25 in jrB| . From (6.4) it follows that 

Lip = -w(x) in D, ip = in D', ip e H? oc (R 3 ). (6.6) 

Thus 

Lip = —w in D, ip = ip N = on S. (6.7) 

Multiply (6.7) by w, integrate over D, then by parts using the boundary conditions (6.7), 
use the equation Lw = and get 



D 



\w(x)\ 2 dx = 0. (6.8) 
Thus w(x) = 0. Estimate (1.18) is proved. □ 
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6.2 Proof of (1.20) and (1.21). 

From (1.19) one gets 

V 2 R + 2i9 • Vi? — q(x)R = q(x) in IR 3 . (6.9) 
Denote L = V 2 + 2i9 • V, and define 



w(x) := L~ L f = -— 3 / 2 ^> - de- (6-10) 



Note that Lw = —f{x). We will prove below that (see |TT], |L5 



||i-VlU^0< c (^) 2 |l/IU 2 W> 9 E M, |0|-oo, (6.11) 

where .Di is an arbitrary compact domain c = c(.Di, ||q , ||l 2 (_b [1 )) ) D C £? a . 
We will also prove that 

||^VlU W < mWfhHD), \9\^oo, 9 EM. (6.12) 

Let us show that (6.11) implies existence of the special solutions (1.19). 
If (6.11) and (6.12) hold, then (1.20) and (1.21) are easily derived. 
Indeed, rewrite (6.9) as 

R = L~ 1 qR + L- l q. (6.13) 
From (6.11) and (6.13) it follows that 



/]oeI0|\* /log 101 V 

L- l qR\\ L ~ iD) < c f -j^ j \\qR\\ LHD) < c ( -^-L j || g || L2(D) ||i2|| Loo(D) . (6.14) 



Therefore the operator L~ x q : L°°(D) — > L°°(D) has the norm going to zero as |0| — ► oo, 

G M. Thus equation (6.13) is uniquely solvable in L°°(D) if |0| > 1, S G M. Moreover, 
the following estimate holds: 

/In |0|\ ^ 

||-R||l°°(d) < c||L _1 g|| LO o (D) < c ( J IMIl 2 (d). (6.15) 

Estimate (1.20) follows. 

To derive (1.21) from (6.12) one writes 

c c 

\\L~ qR\\tf(D) < t^URWlHd) < TZ\h\\w(D)\\R\\L«>(Dy (6-16) 
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Therefore (6.13), (6.15) and (6.16) yield (1.21): 



. |. o || || 2 

\R\\L2(D) < \m\LHD) 



log 






e 





+ \\L 1 q\\L' 2 (D) < 



\9\ 



(6.17) 



Proof of (6.11). \{9 e M then 9 = a + ib, a, b e R 3 , a ■ b = 0, a 2 - b 2 = 1. Choose 
the coordinate system such that a = re 2 , b = te±, r = (1 + 1 2 )^ , ej, 1 < j ' < 3, are the 
orthonormal basis vectors. Then 



e + 29- f = £ + £2 + ^ + 2r6 + 2it& = £ + (6 + r) 2 + £f - r 2 
This function vanishes if and only if 

6 = 0, (6 + r) 2 + e 3 = r 2 . 



2it£ 1 . (6.18) 
(6.19) 



Equation (6.19) defines a circle C T or radius r in the plane 6 = centered at (0, — r, 0). 
Let T$ be a toroidal neighborhood of C T , where the section of the torus by a plane 
orthogonal to C T is a square with size 25 and the center at C T . 

Denote u(x) := L~ l f ', where L" 1 / is defined in (6.10). One has 



(2ny 



Tg e + 2^-9 



+ 



(2nY 



j 

Jr 



rst, i 2 + 2^-9 



h + h, 



(6.20) 



< c 



pdp 



L 1 



t 5 \e + ^-o\ 

-s^Jo ^Ls T^S + (£? + P 2 - r 2 ) 2 

5 ,>:«2 + 2t5+<5 2 ^ 



(6.21) 



'0 J^-2rS+52 a/4£ 2 ^ + yU 2 

<c(L>)||/|| L2(D) / dfi / 
Jo 



_^=,0<*<I 
V 4 ^i 2 + P 2 2 



Li(D) < 



where p 2 = (£ 2 + t) 2 + £ 2 and we have used the Cauchy inequality H/Hl^R 3 ) 
c(D)\\f\\ L 2( D ) and an elementary inequality £ 2 + 2t5 + 5 2 < 3r<5, which holds if £ 2 < 5 2 
and r > 25. 

Let (3:= 2*£l. Then 

- /^d/? ^ 1 '" 3(t+T)a - ^ - 1 



< — 



7T 



dpp ^ dy?- = ^3(t + r)<5 < ctf, (6.22) 



where we have used the relations j 
|0| — > 00. From (6.21) and (6.22) one gets 

h < 4f\\ LHD) 5. 



1 as t — > 00 and took into account that £ — > 00 if 

(6.23) 
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By c > we denote various constants independent of 5 and t. 
Let us estimate I 2 : 



h < C \\J\W 



R3\ T , ie + 2^1 2 



L 2 (D) 



J, 



(6.24) 



where the Parseval equality was used and by J the integral in (6.24) is denoted. One 
has 



J < 



+ 



/ 



dH 



Let us estimate j\- 



■=jl+j2- 



(6.25) 



. ^ r°° „ r pdp 
r ,, r dfi 



< c I — - — — arctg 



2& 



(6.26) 



Let % — x. Then the integral on the right-hand side of (6.26) can be written as: 



c f°°dx 
3i < T / — 

f h x 



7T / r 2 1 



If i — > oo, then ^ — > i. Let us use the elementary inequalities: 



7T 1 7T 

- -x < arctg-, < x < -; 



(6.27) 



(6.28) 



Then 



1 7T X 

arctg- <---, x^+0. 
x 2 2 



7T 1 7T 1 

- - - < arctgy < - - — , y -> +oo. 
2 y 2 2y 



Thus, with A := jp, one gets 



1 



7T / A 

arctg x 

2 6 V x 



1 / A\ c 
<-[x < — , x -> +oo, 



X 



and 



1 



7T / A 

— — arctg x 

2 V x 



<-, x -> +0. 



(6.29) 



(6.30) 



(6.31) 



(6.32) 
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From (6.3), (6.4) and (6.27) one gets 



c 



1 dx \ c , 5 



2/ 



SO 



ln- 

ji < c — 2* , t -> +oo. (6.33) 



Let us estimate j 2 : 

J2 = / d£l I , r ^PP / rf <p 



0<p<T— 



< c 



/ ^ |L (r+5)2 _ r2 4£ 2 /+ p 2 + 4£?f* + 



:2 + 2 _j_ v 2 



(6.34) 



where p = £ 2 + p 2 — r 2 and v = r 2 — p 2 — £ 2 . 
One has: 

5 J£ C2 , ox^ i X2 ^2 C 2 o^X X2 c2 



crd^fir H + 25t + 5 2 t r 2 -i\ t 2t5-5 2 -C 



32 <- / — — arctg h arctg — arctg 



Uo (1 \2 ° 2^ ° 2Cit ° 

c f 5 fn 5 t 5 \ c 

~t Jo "6" U " arCtg 6 + ^6 " arCtg 2^J := 7 J3 ' 



(6.35) 



where we have used the monotonicity of arctg x, for example, 

+ £ 2 + 25r + 5 2 5 
arctg > arctg — , 

etc., and the relation j — > 1 as £ — > +oo, r > t. 
By (6.28), 

7T 5 £i t 71 El 5 IT 2£i . . 

2- arctg- < y , arctg- <---, arctg— >--—. (6.36) 
From (6.35) and (6.36) one gets: 

j, < £ f |i _ li + *A) = £ fx _ I + 2f) = £ ( ! + »V (6 .37) 



Thus 



J 2 < p * -> +oo. (6.38) 
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From (6.20), (6.23), (6.25), (6.33) and (6.37) one gets: 



\u{x)\ < c||/|| z2(D) 



5 + 



In 



it i 



1 

i 



(6.39) 



Choose 8 = |. Then (6.39) yields 

Ms) | < c||/|| L2(D) 

Estimate (6.11) is proved. 

Let us prove (6.12). 

Let L(£) = £ 2 + 20 • £, d = -zV. 



In |6>| 



|0| -»• oo, G M. 



(6.40) 



□ 



Define £(0 := (£ b1 >o T) *• Then 

£(0 = (|£ 2 + 20 ■ £| 2 + 4|£ + 0| 2 + 36)5 > |/ m 0| + 3. 
In 0, vol. 2, p. 31, it is proved that \\L- l f\\ L 2 {Dl) < millf ] cm Wfh^D)- Therefore 



min^ 



min Ce]R 3 | £(0| 



<-^\\f\\^ D ),D(zD u 9eM, |0|-oo, 

1*1 

where c = c(Di, D) > is a constant and we have used the relation 
d|0| < \Im0\ < \9\, ci>0, if 0GM, |0| — »> oo. 
Estimate (6.41) is identical to (6.12). 



(6.41) 



□ 



6.3 Proof of (2.17). 

If p is defined by (2.4), where u(x, a) solves (1.1) then p solves the equation 

V 2 p + 2i6 ■ Vp - q(x)p = q(x) in M 3 , G M. (6.42) 

Let /i = |0|~\ h^0, p(0 := £ 2 + 2(3 ■ f , (3 = hB, (3 ■ (3 = h 2 , \(3\ = 1, 

N := {£ : p(0 = 0, e e K 3 }, ^ := {£ : N)<h,£e M 3 }, iV£ := M 3 \iV ft , 

p = p 1+ ip 2j p 1 = R e p. Note that dP 1 ^ on N, where rfP x is the differential of P v 
Define 



F h u := u 



(27r)2 J R 3 



u 



(x)e' l ^ xh ' 1 dx. 



(6.43) 
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Then 

F h (-id jU (x)) = ihd(.u(Z) = x~u. (6.44) 

Denote 

||p||o := l|p|U 2 (B a ), ||p|| := IMU 2 (R 3 ), \\p\\a,b= \\p\\mB a \B a ), b>a, (6.45) 

\\g(< hD >)p\\ := \\g(VT+e)m)\l D = -iV. (6.46) 

The following Hardy-type inequality will be useful: 
If /(*) G C\-h, h), /(0) = 0, then 

T r 2 |/(0| 2 A < 4 T |/'(t)| 2 dt, > 0. (6.47) 

J-h J-h 

Let us sketch the basic steps of the proof of (2.17) 
Step L If p e C 2 (B r ) and 

P(hD)p:=(hD) 2 p + 2(3-hDp = -h 2 v, v e L 2 (B r ), (6.48) 
where Lq{B 7 ) is the set of L 2 (B r ) functions with compact support in the ball B r , then 

h\\ < hD > 2 p\\ < c\\P(hD)p\\ Vhe(0,h o ), (6.49) 

where ho > is a sufficiently small number. 

Step 2. Let A 1 be a bounded domain with a smooth boundary and A <Z A 1 , rj & 
C yD (A 1 ), < i] < 1, rj(x) = 1 in A, A is a strictly inner subdomain of 

If 

P(hD)p = in Ai, (6.50) 

then 

/i||£>p||l^) < cllplUa^)- (6-51) 

Step 3. Write (6.42) as 

P(hD)p = -h 2 (qp + q). (6.52) 

Let 

V EC™(B b ), 0<r)(x)<l, 
r]{x) = 1 x E B ai , a < cli <b. 
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Then 

P(r]p) = (Pr] - r]P)p - h 2 r](qp + q), P = P(hD). 
Apply (6.49) to (6.53) and get 

h\\<hD > 2 ( pn )\\< C \\(PT,-T,P)p\\ + 

ch 2 \\q\\ L ™ {Ba )\\p\\mB a ) + ch 2 \\q\\ L 2 {Ba) . 
Since rj = 1 in B a , one gets: 

h\\p\\ a < h\\ < hD > 2 M|| < ch 2 \\p\\ a + ch 2 + c\\(Pr] - nP)p\\. 

So 

WpWa^ch + ch-'WiPrj-^pl 

Since = in B a one gets: 

|| (Pr] - rjP) p\\ = \\p(hD) 2 r] + 2h 2 Dr] ■ Dp + 2hp{3 ■ Dr]\\ 
<c(h 2 + h)\\p\\ aub + ch 2 \\Dp\\ auh . 

Using (6.57), one gets 

h\\Dp\\ ai>b < c||p|| 01 _ e)6+e . 

From (6.58), (6.57) and (6.55) one obtains: 

||p|| a <c(/i+||p|| ai _ £ , 6+£ ). 

Since e > is arbitrarily small, the desired inequality (2.17) follows. 
To complete the proof one has to prove (6.49) and (6.51). 

6.4 Proof of (6.49). 

Write (6.49), using Parseval's equality, as 

^ll(l + le| 2 )p||<c||P(Op||. 
IfteNi then h(l + \Z\ 2 )<c\P(Z)\,so 

h 2 [ (i + \z\ 2 ) 2 \m\ 2 dt<c 2 1 \pm 2 \m\ 2 dz 

Jn' Jn: 

h h 

<c 2 I \P(Om\ 2 d£ = c 2 I \P(hD)p\ 2 dx. 
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If £ e N h , then use the local coordinates in which the set N is defined by the equations: 

t = 0, 6 = 0, t = P 1 (0, (6.62) 

and the £i-axis is along vector p defined by the equation f3 — m + ip. Since dP\ 7^ on 
N, these local coordinates can be defined. 

Put / := Pi(Op(0- Then / = at t = 0, / e C°°(R 3 ) if p(x) has compact support, 
and (6.47) yields: 

[ h \p(0\ 2 dt<4 [ h \fl\ 2 dt. (6.63) 

J-h J-h 

Integrating (6.63) over the remaining variables, one gets: 

/ \m\ 2 d£<cf |Va^i(0p(0)| 2 ^<c/ |V € (Pi(0P(O)| 2 de- (6.64) 
Since iV/, is compact, one has 

h 2 [ (i + \er\m\ 2 ^<ch 2 [ \m\ 2 ^. (6.65) 

Using Parseval's equality, S. Bernstein's inequality for the derivative of entire functions 
of exponential type, and the condition supp p{x) C B r , one gets: 

h 2 [ |V 5 (A(0P(0)| 2 ^= / \x\ 2 \P 1 (hD)p(x)\ 2 dx = r 2 [ \P^{hD)p{x)\ 2 dx 

Jm? Jm.3 J B r (6 66) 

< r 2 [ \P{hD)p{x)\ 2 dx. 



From (6.64)-(6.66) it follows that 

h 2 [ (l + |e| 2 ) 2 |p(0| 2 ^<c/ |P(^)p(x)| 2 rfx. (6.67) 
Inequality (6.49) is proved. □ 
6.5 Proof of (6.51). 

Multiply (6.50) by rfp, take the real part and integrate by parts to get: 



h [ 7]\Vp\ 2 dx = ~ [ (pVp + pVp)Vr]dx + 2Re (ify [ pjprjdx 
Ja 1 2 J Al \ J Al 

J \p\ 2 V 2 rjdx + 2Re ^i/3j J pfprjdx^j , 



h 
2 



(6.68) 
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where pj := -j^- and summation is done over the repeated indices. 



Pj ■ d Xj 
One has 



|V 2 ^| < c, 1/3,1 < 1, \2 Pj p\ < ^\ Pj \ 2 + ^\p\ 2 . (6.69) 
From (6.69) and (6.68) one gets: 

hi r]\Vp\ 2 dx < ch I \p\ 2 dx + - [ r]\X7p\ 2 dx + \ I r]\p\ 2 dx. (6.70) 

J A x J Ax 2 J Al h J Al 



Thus 



h 2 I \Vp\ 2 dx <h 2 I 7]\Vp\ 2 dx <c I \p\ 2 dx. (6.71) 

J A J Ax J Ax 

Inequality (6.51) is proved. □ 
Let us prove that 

\\ip(x,0)- u(x,a)v £ (a)da\\ L 2( D) < e, 6 e M, \0\ -> oo, (6.72) 

implies 

H^elU^s 2 ) > ce~ , k = \Im9\, d = diamD, \6\ — > +oo. (6.73) 
Indeed, (6.72), (1.19) and (1.20) imply: 

11/ u(x,a)is £ (a)da\\ L 2 (D) >\\ip(x,9)\\ L 2 {D) -E>ce^, 6 e M, \9\ > 1. (6.74) 

If (6.73) is false for some e > 0, then there is a sequence 9 n e M, — > oo, such that 

IKIU^s^e - ^ -> 0, n -> cx). (6.75) 
This contradicts (6.74) since (6.75) implies 

|| y -u(a;, a)z/ e (o;)(ia|| L 2( £ )) < c||z/ e (a)|| L 2( 5 2) = o ^e^~ j as n — > oo. (6.76) 

Therefore estimate (6.73) is proved. □ 

6.6 Proof of (2.13). 

One has 

/ u(x, a)v(a)da = e l9 ' x (l + p), 
Js 2 
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where 

p := e~ i6 ' x 

"s 2 



/ u(x, a)v(a)da — 1, 
Js 2 

c 

W\' 



1>(x,e) = e»"{l + R), \\R\\ L2{Bbi) <-, 0eM,|0|»l, 



where b\ > b. 

By (1.18), there exist a v(a) such that 



e~ Kbl 



so that 



and 



e««(l + p)- e^(l + R)\\m Bbl) < , « = |/m0|. 

Therefore 

-K&l 

||(p-i?)e^|| L2(Bbi) < 
e- K6l ||p-it!|U( B;)i )< 



ft 
-k6i 



\\p-R\W { b h) <-. 



This implies 



||p||L2(B fcl ) < IIP - R\\L H B bl ) + \\R\\L^B H ) < 7^7 



Thus, inequality (2.13) follows. We claim that ||p||l 2 (B(,) is of order 0(-^). 
Using the above inequalities, one gets: 



p-K&l 

e-^\\p-R\\ LHBh) < \\(p - R)e ie -*\\ L2{Bb) < \\(p - R)e^\\ L , {Bbi) < 



ft 

Thus 

6 -(6i-6)k 
\\p-R\\LHB b )< — ■ 

ri 

Recall that a\9\ < ft < |0|, < c x < \, as \0\ -> oo, 9 e M. Therefore, 

c e _7K 

||p||L 2 (B b ) > PIU W " ||P " R\\L H B b ) — TTjV ~ i 1 = b 1 -b>0. 

Thus, the above claim is verified, since, as \6\ — > oo, 6 G M, one has — — > \[2 and 

Uniqueness class for the solution to the equation Lp = 0. 
Lp := (V 2 + 2«0 • V)p = in M 3 , / |p(x)| 2 (l + |x| 2 )^a; < oo, -1 < £ < 0. (6.77) 
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e -7« 



Taking the distributional Fourier transform of (6.77) one gets: 

£(OP=(£ 2 + 2#-OP = 0. (6.78) 

Thus suppp = C T := {£ : £ G M 3 , L(£) = 0}, and C T is the circle (6.19). By theorem 
7.1.27 in [§, vol 1, p.174] one has: 

/ \p\ 2 ds < c lim sup ( — - / \p{x)\ 2 dx ) . (6.79) 

JC T r ^°° \R J\x\<R J 

Using (6.77) we derive for —1 < £ < 0: 



oo > c> / \u\ 2 (l + \x\ 2 ) i dx > 







2 dx 


(1 + 


X 


2ye\ 



-(ItW/hJ^-^/hJ" 1 ^ 



(6.80) 



Combining (6.79) and (6.80) one gets 



/ \p\ 2 ds < c lim sup ^ = 0, \£\ < 1. 
Jc T # 2 



Thus p(£) = 0, as claimed. □ 
The above argument is valid in M n , n > 2. It was used in El and |TT . 



6.7 Proof of (2.23) 

Let ||^ e || := ||z/ e ||^2(52) and m(e,6) := inf \\u\\ where the infimum is taken over all v G 
L 2 (S 2 ) such that (6.72) holds. 
We wish to prove that 

-bn 

m(e,8) < ce c|e|ln|e| as \9\ -> oo, 9 G M, e= , 6 > a, (6.81) 

K 

where 9 G M, \9\ — > oo, k = \Im9\, and c > stands for various constants. 
Let us describe the steps of the proof. 
Step 1. Prove the estimate 



m(e, 9) < ce Kr — ^ n 2 (e), r > b, 9 G M, e > 0, (6.82) 
\ er J 

where 

ln(n(e)) = ln(| lne|)[l + o(l)], £ -> +0. (6.83) 
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The choice of n(e) in (6.83) is justified below (see (6.97)) and estimate (6.82) is proved 
also below. 

Step 2. Minimize the right-hand side of (6.82) with respect to r > b to get 



m(e,9) < c{2n) n{£) n 2 {e). 

The minimizer is r = 

K 

Step 3. Take e = e(9) = re -> oo, in (6.84). Then 



(6.84) 



In n = ln(nb + In k) [1 + o(l)] = (In re) 



1 + 



In re 



, re — > + 00, 



so, for £ 



one has: 



ci/c < n < c%k, re — > +oo, ci > 0. 
From (6.84) and (6.85) one gets: 

m{9) = m(e(9),9) < ce c ^ ln ^, \9\ -> oo, 9 £ M. 



(6.85) 



(6.86) 



(6.87) 



Estimate (6.81) is obtained. 

Proof o/ fO#J. Since u(a;, a) = (I + T x )e ia ' x where I + T\ is a bijection of C(S 6 ) 
onto C(B b ), inequality (6.72) with D = B b is equivalent to 



(6.88) 



||(/ + Tl )-^_ / e * a - x v £ (a)da\\ LHBb) < ce, 
Js 2 

where c = const > does not depend on e and 9, (I + Ti)~ lr ip — (I + T)ip, 



Tip 



B a 4:7i\x-y 



<i(y)^(y)dy. 



We take b > a, therefore the function ip :— ip + Tip, has the maximal values, as \9\ — > oo, 
of the same order of magnitude as the function ip. The function ip solves the equation 



(V 2 + 1) <p = in R 3 . 

Indeed, (V 2 + l)<p = (V 2 + 1)^ — qip = qip — qip = 0, as claimed. 
Therefore on can write: 



(6.89) 



up := p(x,9) = ^4:7ri e p e Y e (a')j e (r), 



r = \x\, a 



x 
\x\ 



(6.90) 



where Yg are defined in (1.26), ji(r) are defined in (1.29), and ipe = pe(9) are some 
coefficients. 
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It is known that 

oo 

e ia- x = Y^Am e Y^a)Y e (a')j £ (r), (6.91) 



l=o 

SO 



„ oo 

/ e ta - x v £ (a)da = Y j 4m i v £e Y e (a')j e (r), (6.92) 



where v el = (v e ,Y t ) L *(s*)- 
Choose 



v e t = for £ < n(e), v e g = for t > n(e), (6.93) 



where n(e) is the same as in (6.83). 
Then (6.88) implies: 



\(p- le ia - x u £ (a)da\\ 2 L2{Bb) 

oo „5 OO ✓ 7 

= £ leW r*\j e (r)\*dr\<p e \*<c £ 4N 2 (^ 

£=n(e)+l J ° e=n(s)+l ^ 



2^ (6.94) 

< C£, 



u £=n(e)+l 

where formula (1.29) was used. 

From (6.92) and formula (6.91) with a = 6 G M, one gets: 

n(s) n(s) I - 

\wf = E i^i 2 < c E E i^)i 2 ^ cn2 ( £ )rr^YP' Vr > °' (6 - 95) 

e=o E=0m=-e \MeW)\ 

where we have used the formula YH=o Ylm=-e = ( n + 1) 2 > we estimated \ipe\ by the 
coefficient \{e %e ' x , Y e ) L 2^\ 2 = l6ir 2 \Y e (9)\ 2 of the main term of </?, that is, the function 
e l9 ' x , we used estimate (1.28), which gives \(fe\ 2 < c ^^p Vr > 0, and we replaced \je(r)\ 
by \jn( £ ){r)\, the smaller quantity. 

Choose r > b and use (1.29) to get the inequality: 



oo / , \ a oo /u\ n /h\ 2ri ( £ ) 

E MMi) ^ E «*"(;) ^W;J r ><<> < 6 - 96 > 

£=n(s)+l v_ ' £=n(e)+l V / 

,n(e) 



which implies (6.94). Thus (6.94) holds if e Kr (^)™ £ < c^fe, where c stands for various 
constants. One has min r> be Kr = e™ (^)™, and the minimizer is r = ^. Consider 
therefore the equation e n (— ) n = cy^ and solve it asymptotically for n = n(e) as £ — > 0, 
where k > 1 is arbitrary large but fixed. Taking logarithm, one gets In c — \ In ^ = 
n — n In n + n ln(6/t) . 
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Thus | Lne| = In \ = 2ralnra[l + o(l)], and 

In | In e | = (lnn)(l + o(l)), e -> +0. (6.97) 

Hence, we have justified (6.83). 

From (6.94), (6.96) and (1.29), one gets 

e KT (2n(e)) n( - £ ^ 

\\u e \\ < c n 2 (e) — ) ; /; Vr > 6, « = \ImO\, 9 e M. (6.98) 

II £ II — W ( er )"0) ' I I' V / 

Estimate (6.82) is established. □ 
6.8 Proof of (1.30). 

Let Gj be the Green function corresponding to qj(x), j = 1,2. By Green's formula one 
gets 

G 2 (x,y) - Gi(x,y) =/ p(z)G 1 (x, z)G 2 (z,y)dz, p := qi(x) - q 2 (x) . (6.99) 
Take |y| — > oo, — —a and use (1.7) to get: 

u 2 (x, a) — u\(x, a) — / p(z)G\{x, z)u 2 (z, a)dz. (6.100) 

JB a 

Take \x\ — > oo, — a', use (1.7) and (1.2) and get: 

A 2 (a' , a) — A±(a' , a) — — / u±(z, — a')u 2 (z, ct)p(z)dz. (6.101) 

-In.. 

Since A(a',a) = A(—a, —a'), formula (6.101) is equivalent to (1.30). □ 

7 Construction of the Dirichlet-to-Neumann map from 
the scattering data and vice versa. 

Consider a ball B a D D = supp q(x) and assume that the problem 

[V 2 + 1 - q(x)] w = in B a , w = f on S a := dB a , (7.1) 

is uniquely sovable for any / e H%(S a ), where H e (S a ) is the Sobolev space. 
Then the D — N map is defined as 

A : / -> w N (7.2) 
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where wn is the normal derivative of w on S a , N is the normal to S a pointing into 

IK,-- ^\B a . 

If A is known, then q(x) can be found as follows. 

The special solution (1.19)-(1.22) satisfies the equation: 

1>(x) = e t9 - x - [ G(x- y)q(y)il>(y)dy, (7.3) 

JB a 

where G(x) = e ie - x G (x) and V 2 G(x) + G(x) = -5(x) in M 3 . Thus 

V 2 G + 2i6- VG = -5(x), (7.4) 
so that G (x — y) is the Green function of the operator L, see (6.10), that is 

W-wLpTWI- (7 ' 5) 

The function G(x) can be considered known. 
Since qip = (V 2 + one can write, for x G B' a , 

j G(x - s)q(y)i>(y)dy f G(V 2 + l)^dy 

J B a J B a 

= ! [G(x-s)ij N (s)-G N (x-s)ij(s)]ds= j G(x - s)(A - A )ip(s)ds (7.6) 

J S a J S a 

+ / [G(x-s)Aoip-G N (x-s)il>]ds= / G(x - s)(A - Ao)^>(s)ds. 

J S a J S a 

Here A is A for q(x) = 0, we have used Green's formula and took into account that 
/ [G(x - s)A 4j - G N {x - s)ip] ds = j [G(A + l)(p - y?(A + i)G] dy = 0, 

J S a JB a 

where ip solves problem (7.1) with q(x) = and (p — f on S a . 

From (7.3) and (7.6) taking x — > s e one gets a linear Fredholm- type equation 
for ifj\ Sa : 

^ s ) = e ie-s_[ Q{s- s'){K- K Q )i){s')ds'. (7.7) 

JSa 

If A is known, one can find from (7.7) ip\s a and then find q(x) using the following 
calculation. 
Define 

t(9',9):= [ e- w '-yq(y)^(y,6)dy. (7.8) 

JB a 
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By Green's formula, as in (7.6), one gets 

t(0',0)= / e- l6 '- s (A- A )^(s)ds. (7.9) 

JSa 

From (7.8) one gets, using (1.19), (1.20) and (1.9): 

lim t(6', 0)= I e-^ x q{x)dx := q(£). (7.10) 

\9\ — > OO JB a 

6eM 



Therefore the knowledge of A allows one to recover by formula (7.10), but first one 
has to solve equation (7.7). We leave to the reader to check that the homogeneous equa- 
tion (7.7) has only the trivial solution so that Fredholm-type equation (7.7) is uniquely 
solvable in L 2 (S a ) (see a proof in fl33| ). 

Practically, however, there are essential difficulties: a) the function G(x, y) is not 
known, analytically and it is difficult to solve equation (7.7) by this reason, b) the D — N 
map is not given analytically as well. 

Let us show how to construct A from the scattering amplitude A(a', a) and vice versa. 
If A is given then we have shown how to find q(x) and if q(x) found then A(a',a), the 
scattering amplitude, can be found. 

Conversely, suppose A(a', a) is known. Then the scattering solution can be calculated 
in B' a by formula (1.31). 

Let / e Hz(S a ) be given, g(x,y) be the Green function of the operator — V 2 + q(x) — 1 
in IR 3 which satisfies the radiation condition (1.6), and define 

w(x) = / g{x, s)a{s)ds, (7-H) 

JS a 

such that 

w = f on S a . (7.12) 

Since (V 2 + l)w = in B' a , w = f on S a and w satisfies (1.6), one can find w in B' a 
explicitly: 

w ( x ) = X] ]r7~) Y t( a ')ht( r )i r >a, z = \x\, a' = -, (7.13) 
where fe are the Fourier coefficients of /: 

oo 

f(s) = J2fiYi(a'), seS a . (7.14) 
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Therefore the function 

dw(x) 



w N = lim 

\x\-*a,xEB' a Or 

is known. By the jump formula for single-layer potentials one has ( ||16[ , p. 14) 

w% = Wx + a. (7.15) 

The map A : / — > is constructed as soon as we find cr(s), because is already 
found. 

To find a, consider the asymptotics of w(x) as \x\ — > oo, A = [3. Using (1.7) and 
(7.11), one gets: 



(7.16) 



where we have used (7.13) and the asymptotics hi(r) ~ — as r -> +oo. As we have 
already mentioned, the function u(s, a') is known explicitly (see formula (1.31)), and 
equation (7.16) is uniquely solvable for a(s). Analytical solution of equation (7.16) for 
er(s) can be obtained as a series 

oo 

<7(s) = £V / y / (a'), «' = n- (7 - 17) 



1=0 



Substitute (6.91) with a = —j3 into (1.31), take r = a in (1.31) and a' = -, and 
substitute (1.31) into (7.16). By our choice of the spherical harmonics (1.26) both systems 
{Y^}f=o,i,2,... and {Ye}e=o,i,2,... form orthonormal bases of L 2 (S 2 ). Therefore one gets: 



1 °° — . f 

— S^ A-Ki%{-f3)] t {a)a 2 / YAa')a{aa')doi 

H y^^(-/3)/i £ (a)a 2 / lf(a>(aa')(ia' (7.18) 

4?r ^ 2 



£=0 

oo 



fjYejP) 
,._ u Ha) 



E 



Denote 



o-(aa')Y£ jm (a')da' := <j£ m . (7.19) 

S 2 



Using (1.26) one gets: 



Y t ,m(-P) = y, m (-/3) = (-1)%^) 

= (-l) /+/+m *W/?) = (-l) m F,,_ m (/?). 



(7.20) 
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Also define Ai m ^ m > by the formula: 

Ai jm (— /3) = 2_. At >Trh e> m iYii- m t(/3). (7-21) 



I'm! 



The above definition differs from (1.36) and is used for convenience in this section. 
Equating the coefficients in front of Y e _ m (P) in (7.18) one gets 



i e (-l) m j e (a)a 2 a im + — Y] Y] 4m7,m^(a)^m' = rn- ( 7 - 22 ) 

Air / — ' hp(a) 

l>=0m'=-l' ZK ' 

or 

, (-i) m H)^ ^ , . M /^(-l)^-*)' 

**» + 47rj/(fl) W^W^ - a 2 je{a)he{a) ■ ( 7 - 23 ) 

The matrix of the linear system (7.23) is ill-conditioned (see [Q, where estimates of the 
entries of the matrix of (7.19) are obtained and the case of the noisy data is mentioned). 
□ 



Finally let us show (see [33]) that it is impossible to get an estimate 
IK?/||<e(IW||, OeM, H/ll := \\f\\ L * {D) , e(t) - as t - +oo, (7.24) 

if 

Qf= I T(x,y,e)f{y)dy, (7.25) 

J D 

where 

LT := (V 2 + 2i9 ■ V)r = -5(x — y) in D, 6 G M, (7.26) 
r = on S := 9D, (7.27) 

and we assume that 

the problem Lp = 0, p = on S has only the trivial solution . (7.28) 

Indeed, choose a q(x) G L°°(D) such that the problem 

[V 2 + 1 - q{x)] w = in D, w = on S, (7.29) 

has a non-trival solution. 

Define p = e~ %6 ' x w. Then p ^ 0, and 

Lp - qp = in D, p = on 5. (7.30) 

Because of our assumption (7.24), one gets: 

P = / F(x,y)q(y)p(y)dy:=Tp. (7.31) 
Jo 

Were (7.24) true, it would imply for \9\ > 1, 9 G M, that the operator T : L 2 (D) -> 

L 2 (D) in (7.31) has small norm, so p = 0, contrary to our assumption. □ 
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8 Examples of nonuniqueness for an inverse problem 
of geophysics. 

8.1 Statement of the problem. 

In this section the result from |28| is presented. 

Let D C M" := {x : x G W 1 , x n > 0} be a bounded domain, part S of the boundary 
r of D is on the plane x n = 0, f(x,t) is a source of the wavefield, c(x) > is a velocity 



profile. The wavefield, e.g., the acoustic pressure, solves the problem: 

c- 2 {x)u tt - Au = f{x,t) in Dx[0,oo), f(x,t)£0, (8.1) 

u N = on T (8.2) 

M = Mt = at t = 0. (8.3) 



Here iV is the unit outer normal to T, is the normal derivative of u on T. If c 2 (x) 
is known, then the direct problem (8.1)-(8.3) is uniquely solvable. The inverse problem 
(IP) we are interested in is the following one: 

(IP) Given the data u(x,t) Vx 6 S,Vt > 0, can one recover c 2 (x) uniquely? 
The basic result is: the answer to the above question is no. 

An analytical construction is presented of two constant velocities Cj > 0, j = 1,2, 
q / c 2 , which can be chosen arbitrary, and a source, which is constructed after Cj > 
are chosen, such that the solutions to problems (8.1)-(8.3) with c 2 (x) = cj, j = 1,2, 
produce the same surface data on S for all times: 

Ul (x,t) = u 2 (x,t) \/x E S, Vt>0. (8.4) 

The domain D we use is a box: D = {x : aj < Xj < bj, 1 < j < n}. 

This construction is given in the next section. At the end of section 8.2 the data on 
5* are suggested, which allow one to uniquely determine c 2 (x). 



8.2 Example of nonuniqueness of the solution to IP. 

Our construction is valid for any n > 2. For simplicity we take n = 2, D = {x : < 
Xi < 7T, < X2 < 7r}. Let c 2 (x) = c 2 = const > 0. The solution to (8.1)-(8.3) with 
c 2 (x) = c 2 = const can be found analytically 

oo 

u(x,t) = 22 u m(t)(t>m(x), m =(m 1 ,m 2 ) (8.5) 

m=0 
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where 



The data are 



<t>m{x) = -f mim2 cos(mixi) cos(m 2 x 2 ), 
<j) 2 m (x)dx = 1, A0 m + \ m (j) m = 0, 



D 



<pmN = on T, A m := m^ + m 2 , 

1 _ 

700 — — 5 7mi0 — 70m 2 — j 

7T 7T 

7mim 2 = 2/vr if mi > and m 2 > 0, 
c /"* 

u m (t) := u m (t, c) = -= / sin[cV A m (t - r)]/ m (r)dr, 

V JO 



fm(t) := / f(x,t)<fr m (x)dx. 



D 



(8.6) 



(8.6') 



(xi, 0,t) = ^ w m (t, c)7 mim2 cos(mixi). (8.7) 

m=0 

For these data to be the same for c = c\ and c = c 2 , it is necessary and sufficient that 

oo oo 
^2 lm im2 U m (t, Ci) = 7mim 2 Mm(i,C 2 ), Vt > 0, Vmi. (8.8) 

m 2 =0 m 2 =0 

Taking Laplace transform of (8) and using (6') one gets an equation, equivalent to (8.8), 



^ 1 7mim2 f m{p) 



1712=0 



p 2 + c(\ m p 2 + 4\ r 



= 0, Wp > 0, Vmi. 



(8.9) 



Take c\ ^ c 2 , ci,c 2 > 0, arbitrary and find f m (p) for which (8.9) holds. This can be 
done by infinitely many ways. Since (8.9) is equivalent to (8.8), the desired example of 
nonuniqueness of the solution to IP is constructed. 

Let us give a specific choice: c\ = 1, c 2 = 2, f mim2 = for mi ^ 0, m 2 ^ 1 or m 2 ^ 2, 

7M = ^i, 7oi(p) = - (p+ix^+16) - Then ( 8 - 9 ) holds - Therefore, if 



y/2 

f(x, t) = — [/oi(*) cos(x 2 ) + f 02 (t) cos(2a: 2 )] , c ± = 1, c 2 = 2, 

7T 



(8.10) 



then the data Ui(x,t) = u 2 (x,t) Vrr E S,\/t > 0. In (8.10) the values of the coefficients 
are 



/oi(f) = — exp(-i) - ^ 



cos(4t)^sin(4t) 



, f 02 (t) = exp(-t). 



.11) 
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Remark 8.1. The above example brings out the question: 

What data on S are sufficient for the unique identifiability of c 2 (x)? 
The answer to this question one can find in ]76|/ and |71|/ . 

In particular, if one takes f(x,t) = 5(t)5(x — y), and allows x and y run through S, 
then the data u(x,y,t) Wx,y G S, Vi > 0, determine c 2 (x) uniquely. In fact, the low 
frequency surface data u(x,y,k), Vx,y G S Wk G (0, k ), where k > is an arbitrary 
small fixed number, determine c 2 (x) uniquely under mild assumptions on D and c 2 (x). 
By u(x,y, k) is meant the Fourier transform of u(x,y,t) with respect to t. 

Remark 8.2. One can check that the non-uniqueness example with constant velocities 
is not possible to construct, as was done above, if the sources are concentrated on S , that 
is, iff(xi,x 2 ,t) = 5(x 2 )fi(x 1 ,t). 

9 A uniqueness theorem for inverse boundary value 
problem for parabolic equations 

Consider the problem: 



u t + Lu = 0, xED, te[0,T], (9.1) 
u = at t = (9.2) 

u = f(s)5(t) on S. (9.3) 

Here 5(t) is the delta-function, D is a bounded domain in K. n ,n > 3, with a smooth 
boundary S, f G H 3 ^ 2 (S), Lu := —div[a(x)gradu] + q(x)u, a(x) and q(x) are real-valued 
functions, q G L 2 (D), < ao < a(x) < ai, where ao and ai are positive constants, and 
a(x) G C 2 (D), where D is the closure of D. Let h(s,t) := o(s)un, where N is the unit 
exterior normal to S. 

The IP (inverse problem) is: given the set of ordered pairs {f(s),h(s,t)} for all 
t G [0, T], find a(x) and q(x). 

We prove that IP has at most one solution by reducing the uniqueness of the solution 
to IP to the Ramm's uniqueness theorem for the solution to elliptic boundary value 
problem [11]. 

This theorem says: 

Let 

Lu + \u = in D, u = f(s) on S, (9.4) 

and assume that the above problem is uniquely solvable for two distinct real values of 
A. Suppose that the set of ordered pairs {/, h} is known at these values of A for all 
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/ £ H 3 / 2 (S), where h := a(s)u^, and is the normal derivative on S of the solution to 
(9.5). Then the operator L is uniquely determined, that is, the functions a(x) and q(x) 
are uniquely determined. 

We apply this theorem as follows. 

First, we claim that the data h(s, t), known for t £ [0, T] are uniquely determined for 
all t > 0. If 6(t) is replaced by a function rj(t) £ C^°(0,T), f r)(t)dt = 1, then the data 
h(s,t) known for t £ [0,T] are uniquely determined for t > T. 

Secondly, if this claim is established, then Laplace-transform problem (9.1)-(9.3) to 
get the elliptic problem studied in [11]: 

Lv + \v = in D, u = f(s) on S, (9.5) 

and the data H(s, A), where v := / °° e~ xt u(x, t)dt. 

The data H{s, A) := J °° e~ xt h(s,t)dt are known for all A > 0. 

Thus, Ramm's theorem yields uniqueness of the determination of L, and the proof is 
completed. 

We now sketch the proof of the claim: 

The solution to the time-dependent problem can be written as: 

oo 

u(x, t) = Yl e ~ Aj ' c i < ^i( :r )' ( 9 - 6 ) 

where L<pj(x) = Xj(f>j(x) in D,(j)j(x) = on S, ||0j(^)||i,2(d) = 1- The coefficients Cj : = 
- f s f(s)a(s)(f) jN (s)ds. 

Note that the series for u(x,t) and the series obtained by termwise differentiation of 
it with respect to t converge absolutely and uniformly inDx (0, oo), each of the terms 
is analytic with respect to t in the region 5ft > 0, and consequently so are these series. 

Therefore the functions u(x,t) and h(s,t) := UN(s,t) are analytic with respect to t 
in the region > 0, so the data are uniquely determined for t > T as claimed. □ 

At t = the series (9.6) is singular: it does not converge uniformly or even in L 2 (D). 
By this reason the above argument is formal. One can make it rigorous if one replaces the 
delta-function in (9.3) by a C^°(0, T) function r](t), J Q rj(t)dt = 1, and uses the argument 
similar to the one in f52| . 
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